NAVAL 

POSTGRADUATE 

SCHOOL 

MONTEREY,  CALIFORNIA 


THESIS 


FINE-SCALE  VARIABILITY  IN  TEMPERATURE, 
SALINITY,  AND  pH  IN  THE  UPPER-OCEAN  AND  THE 
EFFECTS  ON  ACOUSTIC  TRANSMISSION  LOSS  IN  THE 
WESTERN  ARCTIC  OCEAN 

by 

Stephen  M.  D.  Col 
March  2010 

Thesis  Advisor:  Tim  Stanton 

Co-Advisor:  Daphne  Kapolka 


Approved  for  public  release;  distribution  is  unlimited 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instruction, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send 
comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to 
Washington  headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA 
22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503. _ 

2.  REPORT  DATE 

March  2010 


6.  AUTHOR(S)  Stephen  M.  D.  Col 


11.  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy 
or  position  of  the  Department  of  Defense  or  the  U.S.  Government. 


13.  ABSTRACT  (maximum  200  words) 

Acoustic  propagation  in  the  ocean  is  dependent  on  the  sound  speed  profile,  which  is  a  function  of  temperature, 
salinity,  and  pressure.  In  the  upper  ocean,  these  variables  can  change  quickly  as  a  result  of  fronts,  internal  waves, 
diurnal  heating  and  other  mechanisms.  Transmission  loss  (TL)  is  a  function  of  geometrical  spreading,  absorption,  and 
scattering.  Absorption  is  a  function  of  temperature,  salinity,  pressure,  and  pH.  pH  variability  is  typically  omitted  from 
TL  calculations.  This  thesis  studied  the  effects  of  variability  in  the  vertical  structure  of  temperature,  salinity  and  pH 
on  acoustic  TL  in  the  Arctic  Ocean.  These  parameters  were  measured  in  the  Beaufort  Sea  in  March  2009  during 
ICEX2009  over  a  period  of  15  days. 

Frequencies  considered  in  this  study  were  100Hz,  1kHz,  and  10kHz.  This  range  of  frequencies  brackets  the 
boric  acid  relaxation  frequency  and  is  of  interest  to  the  U.S.  Navy  for  anti-submarine  warfare.  Expected  absorption 
losses  were  estimated  by  integrating  the  loss  over  individual  ray  paths.  Calculations  were  performed  with  constant 
pH  as  well  as  by  using  the  measured  depth  dependent  pH.  The  resulting  difference  in  total  absorption  losses 
computed  with  constant  versus  variable  pH  was  less  than  2dB  over  a  50km  distance.  The  error  in  transmission  loss 
calculations  due  to  ignoring  pH  variability  is  expected  to  be  very  small  compared  to  the  uncertainties  due  to  ice 
scatter-absorption.  Therefore,  routine  measurements  of  the  depth  dependent  pH  levels  are  not  recommended  for  the 
purpose  of  transmission  loss  calculations  under  Arctic  ice  cover. 


16.  PRICE  CODE 


NSN  7540-0 1  -280-5500  Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  239-18 


20.  LIMITATION  OF 
ABSTRACT 


15.  NUMBER  OF 
PAGES 

107 


14.  SUBJECT  TERMS  Acoustic  Propagation,  Transmission  Loss,  Arctic  Ocean,  Temperature 
Salinity  pH  Variability 


18.  SECURITY 
CLASSIFICATION  OF  THIS 
PAGE 

Unclassified 


19.  SECURITY 
CLASSIFICATION  OF 
ABSTRACT 

Unclassified 


17.  SECURITY 
CLASSIFICATION  OF 
REPORT 

Unclassified 


12b.  DISTRIBUTION  CODE 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

9.  SPONSORING  /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

N/A 


5.  FUNDING  NUMBERS 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


4.  TITLE  AND  SUBTITLE  Fine-Scale  Variability  in  Temperature,  Salinity,  and  pH 
in  the  Upper-Ocean  and  the  Effects  on  Acoustic  Transmission  Loss  in  the  Western 
Arctic  Ocean. 


3.  REPORT  TYPE  AND  DATES  COVERED 

Master’s  Thesis 


1.  AGENCY  USE  ONLY  (Leave  blank) 


1 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


11 


Approved  for  public  release;  distribution  is  unlimited 


FINE-SCALE  VARIABILITY  IN  TEMPERATURE,  SALINITY,  AND  pH  IN  THE 
UPPER-OCEAN  AND  THE  EFFECTS  ON  ACOUSTIC  TRANSMISSION  LOSS 

IN  THE  WESTERN  ARCTIC  OCEAN 


Stephen  M.  D.  Col 
Lieutenant,  United  States  Navy 
B.S.,  University  of  Nevada,  Las  Vegas,  2003 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  ENGINEERING  ACOUSTICS 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
March  2010 


Author:  Stephen  Col 


Approved  by:  Tim  Stanton 

Thesis  Advisor 


Daphne  Kapolka 
Co-Advisor 


Daphne  Kapolka 

Chair,  Engineering  Acoustics  Academic  Committee 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


IV 


ABSTRACT 


Acoustic  propagation  in  the  ocean  is  dependent  on  the  sound  speed  profde,  which 
is  a  function  of  temperature,  salinity,  and  pressure.  In  the  upper  ocean,  these  variables  can 
change  quickly  as  a  result  of  fronts,  internal  waves,  diurnal  heating  and  other 
mechanisms.  Transmission  loss  (TL)  is  a  function  of  geometrical  spreading,  absorption, 
and  scattering.  Absorption  is  a  function  of  temperature,  salinity,  pressure,  and  pH.  pH 
variability  is  typically  omitted  from  TL  calculations.  This  thesis  studied  the  effects  of 
variability  in  the  vertical  structure  of  temperature,  salinity  and  pH  on  acoustic  TL  in  the 
Arctic  Ocean.  These  parameters  were  measured  in  the  Beaufort  Sea  in  March  2009 
during  ICEX2009  over  a  period  of  15  days. 

Frequencies  considered  in  this  study  were  100Hz,  1kHz,  and  10kHz.  This  range 
of  frequencies  brackets  the  boric  acid  relaxation  frequency  and  is  of  interest  to  the  U.S. 
Navy  for  anti-submarine  warfare.  Expected  absorption  losses  were  estimated  by 
integrating  the  loss  over  individual  ray  paths.  Calculations  were  performed  with  constant 
pH  as  well  as  by  using  the  measured  depth  dependent  pH.  The  resulting  difference  in 
total  absorption  losses  computed  with  constant  versus  variable  pH  was  less  than  2dB  over 
a  50km  distance.  The  error  in  transmission  loss  calculations  due  to  ignoring  pH 
variability  is  expected  to  be  very  small  compared  to  the  uncertainties  due  to  ice  scatter- 
absorption.  Therefore,  routine  measurements  of  the  depth  dependent  pH  levels  are  not 
recommended  for  the  purpose  of  transmission  loss  calculations  under  Arctic  ice  cover. 
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thennocline  depth  and  the  halocline  depth . 30 

Table  4.  Absorption  Profile  points  of  interest.  Note  that  the  100Hz  and  1kHz  curves 
have  an  absorption  maximum  above  the  sound  channel  and  a  local 
minimum  in  the  channel,  while  the  10kHz  curve  has  a  maximum 
absorption  coefficient  below  the  channel  and  a  local  minimum  near  the 
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I.  INTRODUCTION 


A.  MOTIVATION 

The  Arctic  Ocean  is  an  environment  that  is  acoustically  difficult  to  model. 
Combined  effects  of  growing  and  shrinking  ice,  areas  of  open-ocean,  and  scarcity  of 
ocean  profile  data  yield  significant  errors  in  acoustic  modeling.  This  research  studies 
fine-scale  variability  in  temperature,  salinity,  and  pH  in  the  upper  ocean,  and  their  effects 
on  acoustic  transmission  loss  in  the  Western  Arctic  Ocean.  The  volume  of  sea  ice  is 
declining  faster  than  model  predictions  (Stroeve  et  al.  2007)  and  as  there  are  more  areas 
pennanently  or  seasonally  free  of  ice,  there  will  be  increased  commercial  and  military 
operations  in  the  Arctic  (Beluga  Group  and  Task  Force  Climate  Change  2009).  Enhanced 
understanding  of  the  operational  environment  is  important  at  the  tactical  and  the  strategic 
levels.  With  regard  to  submarine  operations,  application  of  improved  environmental 
sensors  and  models  could  lead  to  increased  perfonnance  of  undersea  weapons  and  higher 
accuracy  in  submarine  detection  ranges.  Since  the  transmission  and  detection  of  acoustic 
energy  in  the  ocean  is  the  primary  method  for  target  prosecution  and  safe  under-ice 
navigation,  the  United  States  Navy  is  keenly  interested  in  effective  and  efficient  use  of 
sonar. 

Considering  that  the  sound  speed  profile  changes  as  a  result  of  tidal,  daily,  and 
seasonal  effects,  as  well  as  longer-scale  climatological  changes,  care  must  be  given  when 
selecting  a  sound  speed  profile  to  use  in  predictive  acoustic  models.  The  primary  goal  of 
this  work  is  to  investigate  the  short-term  effects  of  upper  ocean  variability,  especially  the 
pH,  on  the  transmission  loss  due  to  absorption  between  current  profiles  and  annual 
means. 

B.  BACKGROUND  AND  RELATED  WORKS 

The  undersea  environment  provides  many  challenges  to  navigation, 
communication,  and  combat.  Ocean  water  quickly  attenuates  electromagnetic  radiation, 
including  radar  and  visible  light.  Acoustic  energy,  however,  can  propagate  for  very  long 
ranges  (Kinsler  et  al.  2000).  This  energy  can  be  passively  received  and  processed,  and  in 
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the  realm  of  undersea  warfare  can  be  used  to  detect,  track,  and  classify  vessels.  Sound 
can  also  be  used  actively  and  transmitted  into  the  water  for  communication,  to  localize  an 
enemy  vessel,  or  it  can  be  reflected  from  the  bottom  (or  in  the  Arctic  from  the  ice 
canopy)  to  allow  for  safe  submarine  navigation.  If  a  given  environment  and  the  factors 
affecting  how  sound  travels  in  that  environment  are  known,  detection  ranges  may  be 
more  accurately  predicted. 

The  Arctic  Ocean  is  a  unique  environment  for  submarine  operations.  A  submarine 
must  avoid  time-varying  physical  obstructions  like  icebergs  and  ice  keels  and  deal  with 
an  acoustically  noisy  environment  due  to  cracking  and  shifting  ice.  It  cannot  employ 
radio  communications,  and  it  is  deprived  of  the  classic  submarine  safety  net  of  being  able 
to  surface  and  ventilate  in  an  emergency. 

As  the  Arctic  sea  ice  extent  and  volume  decline  (Stroeve  et  al.  2007),  the  level  of 
interest  in  and  access  to  the  region  increases.  In  September  2009,  two  commercial  ships 
completed  the  first  voyage  from  the  Pacific  Ocean  to  Europe  via  the  Northeast  Passage 
(or  the  Northern  Sea  Route).  The  voyage  from  South  Korea  to  Rotterdam,  Netherlands 
was  about  4,000  km  shorter  than  the  conventional  route  through  the  Indian  Ocean,  Suez 
Canal,  and  the  Mediterranean.  There  is  increasing  interest  in  extracting  the  vast  untapped 
natural  resources  in  the  arctic  seabed  and,  consequently,  effort  is  going  into  resolving 
conflicting  claims  over  Arctic  territory.  As  the  amount  of  traffic  and  commercial  activity 
in  the  region  increases,  there  will  be  an  increased  military  presence  as  well. 

Declining  sea  ice  volume  (Rothrock  et  al.  1999  and  Stroeve  et  al.  2007)  provides 
additional  feedback  effects  for  the  Arctic  ocean-air  interactions.  Thinner  ice  forms  leads 
(breaks  apart)  more  easily.  Reduced  ice  cover  allows  more  summertime  solar  radiation  to 
be  absorbed  by  the  water,  thus  warming  the  ocean  mixed  layer,  and  making  it  more 
difficult  for  ice  to  refreeze  in  the  following  fall  season.  The  resulting  increased  area  of 
open  water  significantly  increases  air-ocean  interactions.  As  atmospheric  carbon  dioxide 
(CO2)  levels  rise  (Bindoff  et  al.  2007  and  Denman  et  al.  2007),  there  is  a  greater  CO2 
transfer  into  the  ocean  (Bates  et  al.2006  and  Bindoff  et  al.  2007).  As  shown  in 
equation(l.l),  when  C02  dissolves  in  water,  it  forms  carbonic  acid,  increasing  the 

concentration  of  hydrogen  ions,  which  lowers  the  pH  of  the  ocean  (Hester  et  al.  2008). 
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The  arrows  pointing  to  the  left  and  right  show  that  the  processes  in  this  weak  acid  system 
can  occur  to  the  left  and  the  right,  depending  on  relative  concentrations  and  external 
stresses  (e.g.,  temperature  or  pressure). 

C02  +  H20  ^  H2CO 3  ^  H+  +  HCO 3  ^  2 H+  +  CO]  (1.1) 

This  lower  pH  reduces  the  amount  of  sound  absorbed  in  the  ocean  (Fisher  and 
Simmons  1977  and  Francois  and  Garrison  1982b),  as  shown  in  Figure  1.  This  will  be 
discussed  further  in  Chapters  II  and  III. 
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Figure  1.  Absorption  vs.  Frequency  (lOFIz  -  lMFIz).  T  and  S  (and  pH  =  7.88)  are 
the  means  from  the  Applied  Physics  Laboratory  Ice  Station  during  Ice 
Exercise  2009  at  10m  depth.  Absorption  values  due  to  reduction  of  pH  0.3  to 

0.5  are  shown. 

1.  Sound  Speed  Profile 

The  deep  ocean  environment  is  typically  slowly  changing  with  range,  and  so  can 
be  considered  range-independent  or  horizontally  stratified  (Lurton  2002).  Acoustic 
energy  travels  as  a  compressional  wave.  If  the  wave  amplitude  and  the  sound  speed  do 
not  change  significantly  over  distances  comparable  to  a  wavelength,  then  the  propagation 
of  acoustic  energy  can  be  modeled  as  following  a  ray  normal  to  the  wave  front  (Kinsler  et 
al.  2000),  with  the  ray  bending  toward  areas  of  lower  sound  speed.  The  speed  of  sound 
traveling  in  the  ocean  depends  on  the  temperature  (T),  salinity  (S),  and  pressure  (P). 
While  there  are  several  methods  for  calculating  sound  speed,  the  Chen-Millero  equation 
(1977),  as  modified  by  Millero  and  Li  (1994),  is  the  standard  for  the  Navy  (Moskal 
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2010).  It  is  referred  to  as  CML  1.1  (for  Version  1.1  of  Chen-Millero-Li).  Though  latitude 
has  been  accounted  for  previously  (Mackenzie  1960),  an  equation  was  recently 
introduced  that  includes  the  correction  factor  for  latitude  and  is  compact  with  a  high  level 
of  accuracy  (Leroy  et  al.  2008). 

2.  Transmission  Loss 

Sound  transmitted  in  the  ocean  spreads  geometrically,  leading  to  lower  power  per 
unit  area.  In  an  isospeed  environment  without  boundaries,  the  geometric,  or  spreading, 
loss  follows  an  inverse  square  relationship  with  distance  from  the  source  (Lurton  2002). 
As  a  sound  wave  travels  in  the  water,  transmission  loss  also  occurs  from  a  variety  of 
mechanisms  including  scattering  and  absorption.  Absorption  results  in  the  loss  of  sound 
energy.  It  is  attributable  to  several  factors  that  include  viscosity,  thermal  conduction, 
thennal  and  chemical  relaxation,  and  interactions  with  inhomogeneities  and  boundaries 
(Kinsler  et  al.  2000).  Scattering  also  results  from  interactions  with  inhomogeneities  and 
boundaries.  However,  scattered  sound  energy  is  not  lost,  it  is  simply  redirected.  These 
mechanisms  will  be  discussed  in  detail  in  Chapter  III. 

C.  APLIS  2009 

The  Applied  Physics  Laboratory  Ice  Station  (APLIS)  was  established  on  the  ice 
cover  in  the  Beaufort  Sea,  350km  north  of  Alaska  in  March  2009  in  support  of  the  U.S. 
Navy’s  Ice  Exercise  (ICEX)  2009.  Personnel  from  the  Arctic  Submarine  Laboratory  and 
University  of  Washington’s  Applied  Physics  Laboratory  (APL-UW)  established  an 
acoustic  tracking  range  and  a  support  camp.  Experiments  were  performed  with  LT  Russ 
Ingersoll,  also  a  student  at  the  Naval  Postgraduate  School  (NPS).  A  working  hut  with  a 
hole  in  the  floor  through  the  ice  allowed  access  to  the  ocean.  Measurements  of  ocean 
temperature  (T),  salinity  (S),  and  pH  were  obtained  with  a  CTD  (conductivity, 
temperature  and  depth  profiler)  and  add-on  pH  sensor.  202  profiles  were  obtained  over 
the  course  of  the  camp.  This  spanned  17  days  while  the  ice  camp  drifted  over  100km.  A 
36-hour  time  series  was  conducted  at  the  beginning  of  the  camp  with  profiles  obtained 
down  to  234m  every  16  minutes.  On  most  days,  six  profiles  were  obtained  to  233m,  but 
over  two  days  four  profiles  were  to  5 1 8m. 
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In  order  to  study  spatial  variability  of  T,  S,  and  pH,  a  radial  survey  was 
conducted,  supported  by  helicopter.  The  helicopter  flew  out  from  the  camp  following  a 
line  of  bearing  of  039°,  and  stopping  every  9.4km,  where  a  hole  was  augured  through  the 
ice  and  two  profiles  were  obtained.  As  a  result,  six  locations  were  sampled,  including  one 
at  APLIS.  After  the  final  radial  was  completed,  47km  to  the  northeast,  a  final  set  of  CTD 
casts  (essentially  a  0th  radial)  was  obtained  from  the  working  hut  at  APLIS.  This  provided 
a  comparison  of  upper  ocean  vertical  structure  with  the  greatest  spatial  difference,  but  a 
small  temporal  difference. 

Preliminary  processing  of  the  data  was  conducted  at  APLIS,  for  diagnostic  purposes. 
This  included  retrieving  the  data  from  the  CTD  and  generating  plots  of  the  temperature, 
salinity,  pH,  density,  and  sound  speed  profiles.  Each  day,  the  sound  speed  profile  was 
provided  to  the  APLIS  Command  Hut,  which  then  relayed  information  to  the  submarines 
involved  in  the  exercise. 

Detailed  processing  and  analysis  was  completed  at  NPS.  During  the  camp,  the  NPS 
students  assisted  a  team  from  the  U.S.  Naval  Academy  and  deployed  two  meteorological 
(ARGOS)  buoys.  Toward  the  end  of  the  camp,  the  students  also  deployed  an  Autonomous 
Ocean  Flux  Buoy  in  support  of  Professor  Tim  Stanton’s  Autonomous  Ocean  Turbulent  Flux 
Buoy  Program  at  NPS,  which  is  supported  by  the  National  Science  Foundation  Office  of 
Polar  Programs. 

D.  ORGANIZATION  OF  THESIS 

This  thesis  will  introduce  how  sound  travels  in  the  Arctic  Ocean,  discussing 
factors  influencing  the  sound  speed  profile  and  factors  influencing  the  attenuation  of 
acoustical  energy.  There  is  much  interest  in  the  Arctic  waters  and  in  the  complex  air- 
ocean  interactions  of  the  region.  Chapter  II  will  review  current  background  research  on 
Arctic  oceanography,  and  the  theory  underlying  acoustic  propagation  and  transmission 
loss  in  the  ocean  will  be  discussed  in  Chapter  III.  Chapter  IV  will  detail  the  experimental 
set  up  and  data  collection  method.  In  Chapter  V  profiles  of  temperature,  salinity,  and  pH 
of  the  Arctic  Ocean  measured  during  ICEX  2009  and  their  sound  speed  structure  will  be 
presented.  Chapter  VI  will  present  conclusions  from  the  research  and  recommendations 
for  application  and  future  study. 
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II.  BACKGROUND 


A.  THE  ARCTIC 

The  Arctic  Ocean  is  a  unique  operating  environment.  Many  regions  of  deep  water 
have  a  sound  speed  profile  that  is  downward  refracting  near  the  surface  and  upward 
refracting  at  deep  depths,  creating  a  thick  channel  where  sound  can  be  trapped  and 
propagate  long  distances.  In  the  Arctic,  the  sound  speed  profile  is  primarily  upward 
refracting,  and  so  sea  ice  reflections  are  expected  to  be  important.  During  the  fieldwork 
conducted  in  the  Beaufort  Sea  in  March  2009,  there  was  a  thin  (100m  deep)  sound 
channel  centered  near  150m,  which  could  allow  extended  propagation  ranges. 

1.  Sea  Ice 

When  seawater  freezes,  it  becomes  first-year  sea  ice.  Icebergs  are  glacial  ice  that 
have  broken  off  into  the  ocean  and  are  not  considered  sea  ice.  In  the  Arctic  Ocean,  some 
sea  ice  melts  during  the  summer,  but  much  of  it  persists  until  the  next  winter,  when  it  is 
considered  perennial  or  multi-year  ice.  In  the  Southern  Ocean,  around  Antarctica,  almost 
all  of  the  ice  is  “seasonal  ice,”  meaning  that  it  melts  and  refonns  annually  (National 
Snow  and  Ice  Data  Center  2010).  The  only  significant  source  of  ice  formation  in  the 
Arctic  summers  is  the  thickening  of  a  “false-bottom”,  a  layer  of  ice  beneath  an  under-ice 
melt  pond  (Notz  et  al.  2003).  Several  factors  influence  the  thickness  of  the  ice,  including 
the  temperature  and  salinity  of  the  ocean  mixed  layer  seawater,  the  amount  of  snow  cover 
on  the  ice,  the  wind  speed,  and  the  air  temperature. 

Sea  ice  that  is  not  “landfasf  ’  is  moved  in  the  ocean,  driven  primarily  by  wind,  but 
also  by  ocean  currents.  While  this  wind  forcing  applies  a  stress  to  the  top  of  the  ice,  the 
ocean  exerts  a  drag  on  the  underside  of  the  ice.  When  diverging  forces  acting  on  an  ice 
floe  exceed  the  tensile  strength  of  the  ice,  it  fractures  and  an  open  body  of  water  is 
formed,  called  a  lead,  that  can  range  in  size  from  centimeters  to  kilometers  (Fricke  1993). 
A  lead  photographed  from  a  helicopter  is  shown  in  Figure  2.  Winter  leads  are  an 
important  source  of  atmosphere-ocean  coupling  (McPhee  and  Stanton  1996),  as  the  rest 
of  the  ice  cover  largely  prevents  any  other  interactions. 
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Figure  2.  Aerial  view  of  a  lead  (left)  that  is  freezing  over,  viewed  from  the  air. 

Close-up  view  of  a  ridge  sail  (right).  Beaufort  Sea.  March  2009. 

If  two  ice  floes  are  forced  together,  their  interaction  is  dependent  on  several 
factors.  Thin  ice  sheets  may  undergo  ‘rafting,’  where  one  sheet  is  pushed  below  the  other, 
and  the  thickness  of  the  ice  where  this  process  occurs  is  effectively  doubled  (Lemke  et  al. 
2007).  Another  mechanism  of  convergence  occurs  when  the  ice  from  either  floe  breaks 
up  along  the  line  where  they  meet.  As  this  process  continues,  some  blocks  are  forced  up 
and  some  forced  down,  forming  the  sail  and  keel  of  a  pressure  ridge.  There  are  several 
manifestations  of  ridges,  but  the  ratio  of  sail  height  to  keel  depth  ranges  from  1/3  to  8 
with  an  average  of  4.0±0.5.  The  ratio  of  keel  draft  to  ridge  width  averages  1/3  (Diachok 
1976).  Regardless  of  the  shape  and  draft,  they  are  features  that  a  submarine  wants  to 
avoid.  Navy  submarines  utilize  upward-looking  sonar  to  ensure  safe  maneuvering  when 
transiting  under  the  ice.  The  U.S.  Navy’s  Arctic  Submarine  Laboratory  (ASL)  has 
archived  ice  draft  observations  since  the  USS  NAUTILUS  (SSN  571)  first  transited  the 
Arctic  Ocean  in  1958  (Submarine  Force  Museum  2009). 

Using  the  submarine  data,  Rothrock  determined  that  the  mean  ice  draft  in  the 
Beaufort  Sea  was  2.1m  (meters)  from  1958-1976,  and  1.2m  during  the  1990s.  The  mean 
ice  draft  trend  for  measurements  from  1993,  1996,  and  1997  showed  a  decrease  of 
0.14m/yr  in  the  Beaufort  Sea  (Rothrock  et  al.  1999). 

In  addition  to  the  ice  thinning,  the  sea  ice  extent  is  also  declining.  Measured  by 
satellite  and  adding  up  areas  with  at  least  15%  sea  ice  concentration,  the  sea  ice  extent 
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has  decreased  (9.12+1.54)%  per  decade  from  1979-2006,  and  it  has  decreased 
(17.91+5.98)%  per  decade  from  1995-2006  (Stroeve  et  al.  2007).  Recent  levels  from  the 
National  Snow  and  Ice  Data  Center  (2010)  are  shown  in  Figure  3.  Model  predictions  and 
data  for  September  are  shown  in  Figure  4.  Those  for  March  are  shown  in  Figure  6.  March 
generally  has  the  greatest  ice  extent  and  thickness,  being  at  the  end  of  the  Arctic  winter, 
whereas  September  toward  the  end  of  the  Arctic  summer  (melt  season)  has  the  smallest 
annual  ice  extent  and  thickness. 


Arctic  Sea  Ice  Extent 
(Area  of  ocean  with  at  least  1 5%  sea  ice) 
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Figure  3.  The  solid  light  blue  line  shows  2009/2010  daily  sea  ice  extent  up  to  March 
4,  2010.  The  dashed  green  indicates  2006/2007  (the  record  low  for  the  month 
of  January);  and  the  solid  gray  indicates  the  average  extent  from  1979  to  2000. 
The  gray  area  around  the  average  line  shows  the  two  standard  deviation  range 
of  the  data.  The  red  ellipse  highlights  the  much  lower  sea  ice  extent  in  the  fall 
compared  to  the  mean.  (National  Snow  and  Ice  Data  Center  2010) 
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Sea  Ice  Extent  (106  km2) 


6  2 

Figure  4.  Arctic  September  sea  ice  extent  (x  10  km")  from  observations  (thick  red 
line)  and  18  IPCC  AR4  climate  models  together  with  the  multi-model 
ensemble  mean  (solid  black  line)  and  standard  deviation  (dotted  black  line). 
Models  with  more  than  one  ensemble  member  are  indicated  with  an  asterisk. 
Inset  shows  9-year  running  means  (Stroeve  et  al.  2007).  Note  that  this  figure 
represents  data  for  September  ice  extent,  two  months  before  what  is  shown  in 

Figure  3. 
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Figure  5.  The  graph  above  shows  daily  sea  ice  extent  for  summer-fall.  The  solid 
light  blue  line  indicates  2009;  the  dashed  green  indicates  2007  (the  record 
low);  and  the  solid  gray  indicates  the  average  extent  from  1979  to  2000.  The 
gray  area  around  the  average  line  shows  the  two  standard  deviation  range  of 
the  data.  (National  Snow  and  Ice  Data  Center  2010) 
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Arctic 

March  Sea  Ice  Extent:  Observations  and  Model  Runs 


,  Ensemble  Mean  ____  Ensemble  ±  aid.  dev. 


. 


9-year  running  means 


Year 


Arctic  March  sea  ice  extent  (x  106  km")  from  observations  (thick  red  line) 
and  18  IPCC  AR4  climate  models  together  with  the  multi-model  ensemble 
mean  (solid  black  line)  and  standard  deviation  (dotted  black  line).  Models 
with  more  than  one  ensemble  member  are  indicated  with  an  asterisk.  Inset 
shows  9-year  running  means  (Stroeve  et  al.  2007). 


2.  Sea  Ice-Albedo  Feedback 

The  ocean  absorbs  sunlight  more  readily  than  ice  does  while  snow  reflects  light 
better  than  ice  or  open  water  so  consequently  the  measure  of  surface  reflectivity,  or 
albedo,  is  an  important  factor  in  sea  ice  growth  and  decline.  The  albedo  is  the  fraction  of 
incident  solar  radiation  that  is  reflected  by  a  surface  (Hanson  1961).  A  good  reflector  has 
a  high  albedo,  and  vice  versa.  The  albedo  of  open  water,  measured  in  the  Beaufort  and 
Chukchi  Seas  in  summer  1958  was  0.04±0.02.  Areas  that  were  80-90%  ice  covered  had 
an  albedo  of  0.40±0.13  (Hanson  1961).  As  the  sea  ice  extent  decreases,  there  is  less  ice 
to  reflect  solar  radiation  away  from  the  Earth’s  surface,  resulting  in  lower  mean  albedo. 
This  allows  more  sunlight  to  be  absorbed  by  the  earth-atmosphere  system,  contributing  to 
a  temperature  increase  (Curry  et  al.  1995).  This  can  make  it  more  difficult  to  form  sea 
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ice,  and  consequently  lowers  the  regional  albedo.  This  is  a  positive  feedback  effect,  and 
so  the  converse  is  also  true  -  as  temperature  decreases,  more  ice  can  form,  increasing  the 
area  albedo,  and  reflecting  more  solar  radiation. 

3.  Air-Ocean  Interactions 

When  the  partial  pressure  of  CO2  in  the  ocean  and  the  atmosphere  are  the  same, 
the  rate  of  CO?  being  taken  up  by  the  ocean  equals  the  rate  of  the  CO2  coming  out  of 
solution,  and  the  system  is  in  equilibrium.  The  following  equation  shows  atmospheric  and 
aqueous  (ocean)  CO2  at  equilibrium,  where  there  is  a  relatively  equal  exchange  to  the  left 
and  to  the  right. 

C02(atm)^  C02(aq)  (2.1) 

The  surface  waters  of  the  Canada  Basin  are  undersaturated  with  respect  to  CO?  in 
the  atmosphere,  and  therefore  have  a  great  capacity  to  take  up  CO2  from  the  air 
(Semiletov  et  al.  2009).  When  the  partial  pressure  of  C02  (pC02 )  in  the  air  increases 
(denoted  by  the  up-arrow  for  atmospheric  CO2)  and  is  greater  than  that  of  the  water, 

co^feco,,,, 

the  equilibrium  function  stresses  toward  right  (denoted  by  the  longer  right  arrow,  and  the 
shorter  left-arrow), 

co„t^cow, 

resulting  in  the  water  taking  up  and  dissolving  more  CO?  from  the  atmosphere, 

co  T — x  co  T 

'^W2(atm)  1  ^  V-'W2(aq)  1  ’ 

and  eventually,  a  new  equilibrium  condition  is  reached 

C02(aq)  . 

The  Revelle  factor  describes  how  pC02  in  the  ocean  changes  given  a  change  in 
dissolved  inorganic  carbon  (DIC).  The  Revelle  factor  is  directly  proportional  to  the  ratio 
of  DIC  to  alkalinity,  and  is  inversely  propotional  to  the  capacity  for  the  ocean  to  absorb 
anthropogenic  CO2  in  the  atmosphere  (Sabine  et  al.  2004).  The  waters  in  the  western 
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Arctic  Ocean  have  an  unusually  high  capacity  to  absorb  CO2.  This  is  reflected  in  the  very 
low  Revelle  Factor  (~3.5-6.5)  for  the  surface  waters  of  the  Canada  Basin,  while  tropical 
and  subtropical  waters  have  Revelle  factors  generally  considered  low  (8-10),  and 
subpolar  and  temperate  waters  have  high  Revelle  Factors  (1 1-14)  (Bates  et  al.  2006).  The 
ability  for  the  waters  of  the  Canada  Basin  to  take  up  so  much  CO  2  is  influenced  by  the 
high  alkalinity  runoff  from  the  Mackenzie  River  and  from  sea  ice  melt  (Bates  et  al. 
2009). 

Even  when  there  is  a  large  ocean  area  covered  by  sea  ice,  there  can  still  be  an  air- 
ocean  exchange  of  CCF  through  leads  and  fractures  (Bates  et  al.  2006  and  Semiletov  et  al. 
2004).  A  circulation  of  CO2  can  be  facilitated  by  density  differences  near  brine  channels. 
The  water  from  a  brine  channel  has  a  higher  salinity  and  density  than  the  surrounding 
water.  Though  the  residence  time  of  surface  waters  in  the  Canada  Basin  is  estimated  to  be 
about  15  years,  these  channels  provide  a  source  for  C02  to  leave  the  upper  ocean  (Bates 
et  al.  2006)  by  entrainment  in  the  sinking  denser  water. 

4.  Ocean  Acidification 

The  equilibrium  equation  for  carbon  dioxide  in  water  is 

c(x  +  h2o  ^  ii2co,  ^  //'  + 1 ico,  ^  21  r  +  co:  (2.2) 

or 

C02  +  4 H20  ^  H2CO,  +  3 H20  ^  H,0  +  HCO,  +  2 H20  ^  2/7,0'  +  CO:  +  H20 
The  second  equation  specifies  that  there  is  excess  water  in  all  stages  of  the  equation. 
When  CO2  dissolves  in  water,  it  forms  carbonic  acid  ( H2CO, ),  bicarbonate  ( /ICO,  ),  and 
carbonate  ( CIO:  ),  and  releases  hydrogen  ions  and  hydronium  ( H 1  and//, <9  ),  which 

lower  the  water’s  pFI  (Hester  et  al.  2008).  This  lower  pH  reduces  the  amount  of  sound 
absorbed  in  the  ocean  (Fisher  and  Simmons  1977  and  Francois  and  Garrison  1982b). 

In  addition  to  increased  CO2  levels,  fossil  fuel  burning  contributes  increased 
sulfur  and  nitrogen  that  cause  ocean  acidification  as  well  (Hester  et  al.  2008),  but  the 
effect  is  small  (Doney  et  al.  2007).  In  the  next  50-100  years,  pH  is  projected  to  decrease 
by  pH  0.3-0. 5  (Bates  et  al.  2009),  and  in  the  next  300  years  by  pH  0.7  (Caldeira  and 
Wickett  2003).  A  simplified  version  of  the  full  Francois-Garrison  absorption  equation 
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(Francois  and  Garrison  1982b)  by  Ainslie  and  McColm  (1998)  shows  how  pH  factors 
into  the  absorption  of  sound  in  seawater. 


a  =  0.106- 


fj2 


fr+f 


(pH- 8) 
,  0.56 


+  0.52(1  +  -^)  A 


43  35  fi+f 


r  rl  -z  T  z  . 

^ +4.9x10 -4/2e 


(2.3) 


In  equation  (3.2),  a  is  the  absorption  coefficient  (dB/km),  f\  is  the  relaxation 
frequency  of  boric  acid  (0.769kHz  for  T=-0.256  °C  and  S=34.67),  fi  is  the  relaxation 
frequency  for  magnesium  sulfate  (41.372  kHz  for  T=-0.256  °C),  and  z  is  depth  in  km. 


Case  1 

Case  2 

Case  3 

f 

0.1  kHz 

1  kHz 

10  kHz 

0.1  kHz 

1  kHz 

10  kHz 

0.1  kHz 

1  kHz 

10  kHz 

z(m) 

1 

0.00144 

0.04946 

1.02269 

0.00088 

0.03314 

0.99994 

0.00065 

0.02624 

0.99032 

30 

0.00145 

0.05032 

1.02746 

0.00089 

0.03365 

1.00410 

0.00065 

0.02660 

0.99423 

70 

0.00149 

0.05742 

1.16074 

0.00092 

0.03832 

1.13249 

0.00068 

0.03025 

1.12055 

180 

0.00089 

0.03840 

1.21621 

0.00057 

0.02763 

1.20015 

0.00044 

0.02307 

1.19337 

450 

0.00114 

0.05152 

1.17101 

0.00071 

0.03490 

1.14416 

0.00053 

0.02787 

1.13281 

1000 

0.00114 

0.04985 

1.11174 

0.00071 

0.03370 

1.08621 

0.00053 

0.02687 

1.07542 

3500 

0.00112 

0.04660 

0.77979 

0.00069 

0.03034 

0.75425 

0.00050 

0.02347 

0.74346 

Table  1.  Absorption  Coefficients  (dB/km)  calculated  using  the  Francois-Garrison 
equation.  Case  1  uses  the  values  in  Table  2.  Case  2  subtracts  0.3  pH  from  the 
values  in  Table  2.  Case  3  subtracts  0.5  pH  from  the  values  in  Table  2.  These 
three  cases  are  demonstrated  graphically  in  Figure  7. 


z(m) 

T(C) 

S 

pH 

1 

-1.43 

26.33 

7.88 

30 

-1.34 

26.62 

7.89 

70 

-0.89 

30.72 

7.95 

180 

-1.53 

32.99 

7.63 

450 

+0.74 

34.83 

7.86 

1000 

-0.05 

34.91 

7.85 

3500 

-0.28 

34.95 

7.87 

Table  2.  Mean  values  of  T,  S,  and  pH  at  the  given  depths.  From  APLIS  2009. 
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Absorption  (dB/km) 


Figure  7. 
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Absorption  \s.  Depth.  Demonstrate  lowering  absorption  coefficient  with  lowering  pH 


0  500  1000  1500  2000  2500  3000  3500 


Depth  (m) 

Graphical  representation  of  data  in  Table  1.  The  frequency  for  the  three 
sets  of  lines  are  10kHz  for  the  top  set,  1kHz  for  the  middle  set,  and  100Hz  for 
the  bottom  set.  In  each  frequency  set,  the  green  lines  represent  the  measured 
pFI  values  at  APLIS  2009.  The  blue  lines  represent  pFI  lowered  by  0.3  pH. 
The  red  lines  represent  pH  lowered  by  0.5  pH. 
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III.  THEORY 


A.  ACOUSTIC  TRANSMISSION  LOSS  IN  THE  OCEAN 


Acoustic  Transmission  Loss  (TL),  is  defined  e.g.,  (Kinsler  et  al.  2000)  as  shown 
in  equation(3.1).  This  shows  the  relationship  between  P(r),  the  pressure  amplitude  at  a 
range  r  from  a  sound  source,  and  P  (1),  the  pressure  amplitude  at  lm  from  the  source. 
Transmission  loss  terms  are  often  grouped  into  categories  of  absorption,  geometrical 
spreading,  and  scattering  effects,  as  shown  in  equation(3.2). 


TL  =  10  log 


m 

yl(r) 


=  20  log 


Pj 1) 
\P(rb 


tt  -7T  j-  77  +  77  +  TT 

1  ^  1  ^ absorption  ^  1  ^ geo  ^  1  ^ bottom  ^  1  ^ surface 


(3.1) 

(3.2) 


1.  Absorption 


Many  factors  affect  the  absorption  of  acoustic  energy  in  seawater,  including 
thermal  conduction  and  thermal  relaxation  (Kinsler  et  al.  2000),  the  primary  mechanisms 
are  viscous  absorption  (significant  above  100kHz),  and  chemical  relaxation  due  to  ocean 
salts,  mainly  magnesium  sulfate  and  boric  acid  (Ainslie  and  McColm  1998  and  Francois 
and  Garrison  1982b).  The  Francois-Garrison  equation  for  total  absorption  is  given  in 
equation  (3.3).  Above  hundreds  of  kilohertz,  the  pure  water  effect  dominates.  The 
magnesium  sulfate  effect  dominates  in  the  region  from  a  few  kilohertz  to  hundreds  of 
kilohertz.  Below  a  few  kilohertz,  the  primary  effect  is  the  chemical  relaxation  of  boric 
acid  (Kinsler  et  al.  2000). 


AxPJxf2  A2P2fJ2  A  _ 
a=  \ 1  1  ,  +  2~  2  ,  +A,PJZ 

f2+f{  f2+f: 2 


(3.3) 


In  equation(3.3),  the  subscript  ‘1’  are  for  boric  acid,  ‘2’  is  for  magnesium  sulfate,  and  ‘3’ 
is  for  pure  water.  In  this  equation,  ^4v  is  a  coefficient  that,  for  boric  acid  depends  on 


temperature  (T,  °C),  salinity  (S,  %o,),  pFI,  and  depth  (D,  m).  (Salinity  used  to  be  measured 
by  determining  the  mass  (kg)  of  salt  per  kg  of  seawater  (given  as  parts  per  thousand 
(ppt,%o),  assuming  that  the  composition  of  sea  salt  is  constant.  Later  it  was  determined  by 
a  conductivity  ratio  and  listed  as  %o,  but  it  had  problems  with  ambiguities  due  to  ionic 
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ratio  variations.  The  introduction  of  the  Practical  Salinity  Scale  of  1978  calculates 
salinity  as  a  unitless  conductivity  ratio  to  a  standard  KC1  solution  (Siedler  and  Peters 
1986).  The  numerical  difference  between  practical  salinity  and  ppt  is  small  (Sea-Bird 
Electronics  2010).  Equation(3.3)  is  implemented  with  practical  salinity.  For  magnesium 
sulfate,  Ax  depends  on  T,  S,  and  depth;  and  for  pure  water  depends  on  T  only.  P  is  the 
pressure  dependency,  fx  is  the  relaxation  frequency,  and  /  is  the  frequency  of  sound  that 

is  of  interest.  In  equation  (3.3),  the  subscript  ‘1’  is  for  boric  acid,  ‘2’  is  for  magnesium 
sulfate,  and  ‘3’  is  for  pure  water.  The  boric  acid,  magnesium  sulfate,  and  pure  water 
dependencies  are  detailed  in  the  following  sections. 

a.  Pure  Water 

The  amount  of  sound  absorption  in  pure  water  is  a  function  of  the 
temperature,  pressure,  and  the  frequency  of  the  sound.  In  general,  as  the  temperature  goes 
up,  the  absorption  goes  down.  Increasing  depth  has  a  similar  effect,  where  deeper  depths 
have  lower  absorption  values  (Francois  and  Garrison  1982a).  Higher  frequencies  have 
exponentially  higher  absorption  (Francois  and  Garrison  1982a  and  Kinsler  et  al.  2000). 

For  temperatures  less  than  or  equal  to  20  °C,  the  pure  water  contribution  to 
the  absorption  shown  in  (3.3)  is 

A  =4.937x10^-2.59x10  5r  +  9.11xl0'7r2  -1.50x10  8T3  ( — dB—) 

km -kHz2 

P3  =  1 -3.83  x  1CT5Z)  +  4.9  xlCT10!)2.  (3.4) 

b.  Magnesium  Sulfate 

The  sound  speed  at  a  given  point  is  a  function  of  temperature,  salinity,  and 
depth,  and  the  absorption  of  sound  due  to  magnesium  sulfate  ( MgSO 4 )  is  a  function  of  T, 
S,  depth,  and  sound  speed  (Francois  and  Garrison  1982a).  The  contribution  of  a  chemical 
relaxation  that  involves  magnesium  and  carbonic  acid  (Mellen  et  al.  1979)  is  small  and  is 
considered  to  be  included  in  first  two  terms  of  equation  (3.3).  The  equations  (3.5)  show 
the  MgSO 4  contribution  to  the  total  absorption. 
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(3.5) 


=  2 1 .44-  (1  +  0.0257)  ( — — — ) 
c  km  ■  kHz 


R  =1-1.37x10  D  + 6.2x10  D 


9  r,2 


/2  = 


8.17x10 


(8-1990, 


1  +  0.001 8(5  -35) 


{kHz) 


In  equations  (3.5)  and(3.6),  cis  the  speed  of  sound  (m/s),  and  #  =  273  +  7  is  absolute 
temperature  in  Kelvin. 


c.  Boric  Acid 

The  effect  that  the  relaxation  of  boric  acid  has  on  sound  absorption  is 
interesting  in  that  it  is  dependent  on  the  pH  of  the  water  in  addition  to  the  temperature, 
salinity,  depth,  and  sound  speed  (Fisher  and  Simmons  1977  and  Francois  and  Garrison 
1982b).  As  pH  goes  up,  the  absorption  increases,  and  this  relationship  is  exponential 
(Kinsler  et  al.  2000  and  Francois  and  Garrison  1982b).  This  is  shown  in  equation  . 

h,=^xio<o.W  (  dB  ) 
c  km  ■  kHz 

P,=  1  (3-6) 

rrr  1245 

/i=  2.8^10"'  {kHz) 

2.  Geometrical  Spreading 

An  accurate  prediction  of  geometrical  transmission  loss  using  ray  theory  requires 
a  model  that  sums  up  the  contribution  of  sound  energy  from  source  to  receiver  along  all 
possible  paths  by  keeping  track  of  the  spreading  of  each  path’s  ray  bundles. 
Unfortunately,  although  such  models  are  available,  none  was  found  that  could  handle  the 
variable  pH  effect  on  absorption.  Furthermore,  since  the  focus  of  this  work  was  mainly 
on  the  question  of  whether  pH  was  important  to  measure,  the  decision  was  made  to  look 
at  the  effect  of  pH  on  the  absorption  losses  for  each  path  individually  and  to  assume 
spherical  spreading  for  the  geometric  loss.  This  approach  yields  a  rough  guess  of  the 
geometrical  transmission  loss  along  each  path  that  is  independent  of  pH.  Therefore,  the 
effect  of  absorption  due  to  pH  is  preserved  along  these  paths. 
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If  sound  energy  radiates  away  from  a  source  spherically,  the  surface  area  of  the 
wave  front  increases  proportional  to  the  growing  surface  area  of  a  sphere,  as  a  function  of 
the  radius  squared.  It  follows  that  its  pressure  amplitude  decreases  accordingly,  and  so 
the  TL  due  to  spherical  spreading  can  be  written  TLsphencai  =  20  log  r  (Kinsler  et  al.  2000). 

This  thesis  is  focused  on  propagation  ranges  on  the  order  of  50  km.  Due  to  the 
strong  upward  refracting  environment,  sound  rays  that  leave  a  source  at  steep  launch 
angles  (greater  than  17°)  incur  additional  losses  from  interactions  with  the  bottom  and  are 
assumed  to  have  a  negligible  contribution  to  the  signal  at  long  ranges  (Stotts,  S.  A.,  R.  A. 
Koch,  and  N.  R.  Bedford  1994)  (see  also  the  next  section). 

3.  Bottom  Loss 

Depending  on  the  bottom  structure  and  composition,  acoustic  energy  can  be 
reflected,  scattered,  and  absorbed.  The  bathymetry  in  the  Beaufort  Sea  where  APLIS 
2009  was  established  was  very  flat  and  roughly  3500m  deep  (Jakobsson  et  al.  2008).  The 
bottom  type  in  this  area  is  fine  silt  and  mud  (Hunkins  et  al.  1970).  Launch  rays  were 
considered  (+15°)  that  do  not  interact  with  the  bottom,  since  steep  angle  rays  that  interact 
with  the  bottom  undergo  additional  losses  and  will  not  have  significant  influence  on  long- 
range  observations  (Stotts  et  al.  1994),  and  so  TLboUom  =  0 . 

4.  Surface  Loss, 

a.  Sea  Surface  Effects 

The  sea  surface  is  a  pressure-release  boundary,  and  so  sound  which  is 
incident  upon  it  is  reflected  with  a  phase  shift  of  n.  The  roughness  of  the  sea  surface  is  a 
function  of  the  sea  state,  which  is  driven  by  wind.  Wave  action  can  entrain  bubbles, 
which,  for  acoustic  frequencies  above  those  on  the  order  of  10kHz  the  primary  source  of 
backscattering  (Dahl,  P.  H.,  J.  W.  Choi,  and  H.  C.  Graber  2008).  There  also  tends  to  be 
more  biologic  activity  near  the  surface,  which  can  serve  to  scatter  incident  acoustic 
energy.  During  the  course  of  the  data  collection,  March  2009  in  the  Beaufort  Sea,  the 
open-water  fraction  was  very  low,  and  so  the  sea  surface  effects  were  assumed  to  be 
small,  and  dominated  by  the  ice  effects. 
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b.  Ice  Effects 

Acoustic  energy  is  scattered  by  interactions  with  under  ice  features.  Two 
primary  types  of  features  are  ice  edges  and  ridges  (Fricke  1993).  Ice  edges  occur  at  the 
boundaries  of  leads  and  polynyas  (regions  of  no  ice  surrounded  by  ice  coverage  (Roth 
2003)).  If  an  acoustic  wave’s  angle  of  incidence  on  the  ice  is  great  enough,  it  can  get 
transmitted  into  the  ice,  travel  through  the  ice,  and  get  reflected  off  of  the  air-ice 
interface.  Work  in  this  area  was  done  using  geophones  buried  in  the  ice  to  measure 
sound  coupled  from  the  ocean.  It  was  shown  that  for  long  range  propagation,  the  rays  that 
penetrate  the  ice  are  not  important  when  the  source  and  receiver  are  both  in  the  water 
(Stotts  et  al.  1994).  In  order  for  an  acoustic  ray  to  penetrate  into  the  ice,  it  needs  to  have  a 
steep  angle  of  incidence.  This  results  from  the  much  larger  speed  of  sound  in  ice, 
compared  to  the  speed  of  sound  in  the  seawater  below  the  ice.  This  difference  in  sound 
speed  leads  to  a  large  difference  in  the  indices  of  refraction. 

Consider  an  acoustic  ray  that  is  incident  from  the  water  to  the  ice,  where 
the  ice  has  sound  speed  Ci =  3,200  m/s,  (Kinsler  et  al.  2000),  and  the  water  has  Ci  =  1,430 
m/s.  Then  the  index  of  refraction  n=ci/c2.  The  minimum  (critical)  0C  required,  for  a  ray  to 
propagate  in  the  ice  is  0C  =  cos'*(n)  =  63.5°.  Judging  from  this,  all  of  the  rays  within  an 
angle  of  +/- 15°  are  reflected  back  into  the  water  and  do  not  penetrate  the  ice. 


The  scatter-absorption  contribution  to  the  transmission  loss  due  to 
interactions  with  the  ice  is  calculated  from  TLice  =  -lOlogR ,  and  Diachok’s  1976 


equation  (3.7),  which  gives  the  energy  reflectivity  R.  R  is  the  ratio  of  energy  that  gets 
reflected  cleanly  in  the  initial  travel  direction.  The  acoustic  ray  is  assumed  to  be 
otherwise  perfectly  reflected  from  the  ice.  For  this  equation,  ridges  are  assumed  to  be  half 
cylinders,  with  length  much  greater  than  width,  with  average  keel  depth  d,  and  spatial 
frequency  N.  x  is  the  correction  for  the  eccentricity  of  an  ellipsoid  ice  ridge  form.  This  is 
assumed  to  be  1.0  since  Diachok  assumes  a  ridge  is  half  of  a  right  circular  cylinder. 


(3.7) 
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This  theoretical  equation  proposed  for  frequencies  such  that  kd>l,  where  k  is  the 
wavenumber:  k=2*n [/A.  For  the  ice  drafts  and  frequencies  studied,  kd>l,  and  so  this  ice 
scatter-absorption  model  is  independent  of  further  considerations  of  frequency. 

B.  SOUND  SPEED  PROFILE 

Sound  travels  through  a  medium  as  a  compressional  wave  in  a  series  of 
compressions  and  rarefactions.  The  sound  speed  profde  is  often  considered  to  be  radially 
homogeneous  and  range-independent.  Even  so,  inhomogeneities  in  temperature  and 
density  change  how  a  sound  ray  refracts  as  it  propagates.  This  study  considers  vertical 
structure  as  measured  by  CTD  in  the  upper  ocean  and  concatenated  with  deeper  data 
sources  to  develop  sound  speed  profiles  for  the  full  water  column.  It  will  not,  however, 
consider  range-dependence  of  the  sound  speed  profile  because  over  the  50km  range 
surveyed,  the  variability  in  SSP  due  to  spatial  separation  was  not  significantly  different 
from  the  variability  due  to  the  ice  sheet  drifting  or  the  time  separation  between  profiles. 

Similarly,  differing  media  can  have  quite  different  sound  speeds.  For  example,  the 
speed  of  sound  in  air  at  0°C  is  331.5m/s,  in  molecular  hydrogen  at  0°C;  it  is  l,269.5m/s. 
In  steel,  the  speed  is  nominally  6, lOOm/s,  but  only  3,200m/s  in  ice  (Kinsler  et  al.  2000). 
Near  the  surface  in  the  Beaufort  Sea,  the  annual  average  seawater  sound  speed  is 
l,434m/s,  and  l,507.5m/s  near  the  bottom  at  3,500m.  These  sound  speeds  were 
calculated  from  World  Ocean  Atlas  2001  salinity  data  (Boyer  et  al.  2002)  and 
temperature  data  (Stephens  et  al.  2002).  The  speed  of  sound  in  the  ocean  has  been 
measured  extensively.  There  is  little  ambiguity  about  the  measured  sound  speed  for  a  set 
of  given  conditions.  There  is,  however,  much  debate  about  how  to  develop  a 
mathematical  equation  to  describe  how  the  speed  of  sound  in  the  ocean  varies  as  a 
function  of  temperature,  salinity,  and  pressure  (which  is  related  to  depth  and  latitude) 
(Mackenzie  1960,  Wilson  1960,  Del  Grasso  and  Mader  1972,  Del  Grasso  1974,  Chen  and 
Millero  1977,  and  Leroy  et  al.  2008). 

The  Navy’s  Oceanographic  and  Atmospheric  Master  Library  first  accepted 
Wilson’s  1960  equation  as  the  Navy  Standard.  When  the  Navy’s  mission  focus  shifted  to 
more  littoral,  “brown  water”  areas,  many  sound  speed  equations  were  tested.  The 
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approved  standard  sound  speed  algorithm  for  the  United  States  Navy  is  the  Chen  and 
Millero  equations  from  1977,  as  modified  by  Millero  and  Li  in  1994.  Recently,  however, 
there  were  issues.  The  algorithm,  implemented  as  CML  1.0,  would  give  erroneous  SSP 
results  when  users  attempted  to  calculate  sound  speed  profiles  in  regions  like  the  Red 
Sea.  This  was  because  there  were  instances  when  the  salinity  (>40)  was  higher  than  the 
limits  in  the  program.  The  previous  limits  were:  temperature  (0  to  40)  °C,  salinity  (5  to 
40)  ,  and  pressure  (0  to  1,000  bar)  (Krynen  2010).  The  limits  were  expanded  and  the 
modified  algorithm,  CML  1.1,  was  approved  January  26,  2010  as  the  new  Navy  Standard 
Sound  Speed  Algorithm  (Moskal  2010).  The  revised  limits  are:  temperature  (-3  to  40)  °C, 
salinity  (0  to  45)  ,  and  pressure  (0  to  1,100)  bar  (Naval  Oceanographic  Office  2009).  A 
comparison  of  SSPs  generated  by  several  equations  is  shown  in  Figure  8.  Given  the 
profiles  obtained  at  APLIS  2009,  the  difference  between  algorithms  is  insignificant.  The 
maximum  difference  is  at  3500m,  when  there  is  a  0.009%  difference  between  CML  1.1 
and  Mackenzie,  1960,  as  shown  in  Figure  8.  In  the  upper  ocean,  the  variability  in  SSP 
due  to  selection  of  equation  set  is  dominated  by  variability  due  to  unique  T  and  S 
profiles. 
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Figure  8.  Comparison  of  sound  speed  profile  equations.  The  plot  on  the  left  uses  the 
listed  algorithms  for  determining  sound  speed  from  all  the  profiles  of  T  and  S 
at  APLIS  2009.  The  plot  on  the  right  uses  the  same  data  and  is  zoomed  into 
the  region  where  most  data  acquisition  occurred.  Note  that  profiles  were 
conducted  to  a  maximum  depth  of  5 1 8m.  The  variability  due  to  algorithm 
choice  in  the  upper  ocean  is  dominated  by  the  consequences  of  T  and  S 

variability. 
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IV.  EXPERIMENTAL  SET  UP  AND  DATA  COLLECTION 


A.  EQUIPMENT 

1.  Seabird  SBE19plus  V2 


Figure  9.  Sea-Bird  Electronics  SBE  19plus  V2  with  SBE18  pH  sensor.  As  deployed 

at  APLIS  2009. 

This  conductivity-temperature-depth  profiler  (CTD)  employed  a  SBE  5M  Mini 
Submersible  Pump  which  maintained  a  constant  flow  through  the  conductivity  cell  and 
temperature  sensor  regardless  of  CTD  motion.  It  minimized  salt  spiking  arising  from 
mismatched  time  responses  between  the  temperature  sensor  and  conductivity  cell.  The 
resulting  fast  response  along  with  a  slow  vertical  travel  time  of  0.5m/s  provided  enhanced 
resolution  of  the  vertical  structure.  Sampling  at  4Hz,  eight  measurements  were  logged 
every  meter.  In  the  Arctic  environment,  it  was  very  important  to  protect  the  glass 
conductivity  cell  from  damage  and  freezing.  If  water  freezes  onto  a  conductivity  cell,  the 
measurements  are  inaccurate.  It  was  also  important  to  ensure  that  the  water  was  clear 
prior  to  lowering  the  CTD  into  the  water  to  prevent  debris  from  entering  the  instruments. 
When  the  CTD  was  not  in  use,  it  was  suspended  in  the  hut,  where  the  temperature  was 
maintained  near  22°C. 

As  an  add-on  to  the  CTD,  a  SBE  1 8  pH  sensor  was  used  to  study  the  structure  and 
variability  in  the  pH  of  the  upper  ocean.  This  sensor  uses  a  pressure-balanced  glass- 
electrode  Ag/AgCl-reference  pH  probe.  For  most  of  the  profiling,  it  was  attached  to  the 
side  of  the  CTD  housing,  as  shown  at  the  bottom  of  Figure  9.  In  this  configuration,  the 
diameter  of  the  CTD  was  0.23m.  The  largest  hole  that  could  be  drilled  with  the  portable 
auger  was  0.1524m,  so  for  the  helicopter-assisted  radial  casts,  the  sensor  was 
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repositioned  and  attached  to  the  bottom  end  of  the  CTD  housing.  This  repositioning 
lowered  the  pH  sensor  0.5m,  and  the  resulting  CTD  diameter  was  about  0. 15m. 

The  intent  had  been  to  profile  to  300m,  but  upon  testing  the  winch  and  pulley 
apparatus  with  a  7kg  test  weight  (CTD  as  used  was  8.3kg),  it  became  evident  that  instead 
of  305m  of  line  there  was  only  about  250m  of  line.  It  was  decided  to  profile  to  233m  to 
allow  a  safety  buffer  of  line  on  the  spool.  A  Shimano  Dendou-Maru  deep-sea  fishing  reel 
was  used  as  the  winch,  with  a  very  thin  Spectra/Nylon  100  lb-test  line  going  through  a 
small  pulley  system. 

B.  DATA  COLLECTION 

1.  Work  Area  Setup 

Most  of  the  CTD  operations  were  performed  from  a  working/sleeping  hut 
assembled  by  the  APL-UW  APLIS  team.  The  hut  was  constructed  on  1.28m  of  ice,  and 
placed  over  a  0.46m  hydro-hole  that  was  melted  through  the  ice  sheet. 

2.  Equilibration 

Initially,  the  CTD  was  allowed  to  equilibrate  for  eight  minutes  at  lm  but  it  was 
later  determined  that  it  could  be  more  effective  to  lower  the  CTD  to  7.5m  for  15  minutes 
for  equilibration.  The  reason  for  lowering  the  CTD  further  for  equilibration  was  that  the 
hydro-hole  was  maintained  open  by  two  heat  tapes  on  opposite  sides  of  the  hole, 
extending  to  the  bottom  of  the  ice.  In  addition,  the  surface  was  kept  from  freezing  by  a 
fan  circulating  warm  air  from  the  top  of  the  hut  to  the  surface  of  the  hole.  Lowering  the 
CTD  to  7.5m  allowed  the  sensors  to  equilibrate  away  from  the  warmer,  fresher  water  in 
the  hole.  It  was  subsequently  noted  that  there  was  little  change  in  the  parameters  after 
five  minutes,  and  the  final  equilibration  routine  was  established  as  7.5m  for  five  minutes. 

3.  Time  Series 

To  study  the  tidal  and  diurnal  variability  in  T,  S,  and  pH,  a  36-hour  time  series 
was  conducted.  During  the  time  series,  124  profiles  were  obtained  to  an  average  depth  of 
233m,  lowering  the  CTD  every  14.5  minutes.  While  continuous  profiling  for  the  entire 
time  would  have  been  desirable,  the  pulley  wheel  often  seized  and  required  periodic 
maintenance  and  lubrication. 
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4.  Daily  Casts 

After  the  time-series,  six  CTD  profiles  to  233m  were  obtained  each  day,  one 
profile  or  down-cast  for  the  descent,  and  one  profile  or  up-cast  for  the  ascent.  The  sound 
speed  profiles  were  calculated  and  transmitted  to  the  participating  submarines  each  day 
so  that  they  could  more  accurately  predict  acoustic  propagation  in  the  ocean.  Twice  the 
water  was  profiled  to  5 1 8m  making  use  of  a  heavy-duty  winch  assembly  borrowed  from 
APL-UW. 


5.  Radial  Survey 

In  order  to  study  the  spatial  variability  of  the  water  properties  and  therefore  the 
spatial  variability  and  range  dependence  of  the  sound  speed  profile,  a  survey  away  from 
the  ice  camp  was  conducted.  Using  a  helicopter  for  transportation,  two  profiles  were 
obtained  at  each  of  five  locations.  The  average  spacing  was  9.4km,  and  the  line  of 
bearing  from  the  camp  was  039°.  Including  two  profiles  conducted  back  at  the  camp  (the 
0th  radial  stop),  twelve  profiles  were  obtained  in  all. 

Three  of  the  down-casts  from  field  surveys  were  erroneous  due  to  a  suspected  ice 
clog  on  the  intake  of  the  conductivity  cell.  The  temperature  and  salinity  readings  returned 
to  normal  halfway  down,  but  because  of  this  problem,  these  three  profiles  were 
subsequently  discarded  and  replaced  with  the  corresponding  up-casts.  The  up-casts  and 
down-casts  were  very  similar. 

C.  DATA  UPLOAD/CONVERSION 

After  uploading  the  data  from  the  CTD,  it  was  converted  from  .hex  format  with 
SBE  Data  Processing,  Version  7.17a  (Sea-Bird  Electronics  2008).  The  temperature  (T)  is 
output  as  degrees  Celsius  according  to  the  International  Temperature  Scale  of  1990  (ITS- 
90).  The  conductivity  is  output  in  mS/cm  (milliSiemen/centimeter).  The  pressure  is 
converted  and  output  as  depth  in  meters  of  seawater.  pH  is  calculated  from  the  voltage  of 
the  pH  sensor  and  the  temperature.  Salinity  was  then  calculated  by  the  equations  from  the 
Practical  Salinity  Scale  of  1978  (PSS-78)  and  is  expressed  as  a  unitless  conductivity  ratio 
(Fofonoff  and  Millard  1983). 
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D.  EQUIPMENT  MAINTENANCE 

After  profiling,  the  entire  CTD  was  rinsed  with  fresh  water.  The  conductivity  cell 
was  then  flushed  with  0. 1-0.2%  Triton  X-100  solution,  deionized  (DI)  water  or  fresh 
water,  and  drained.  Prior  to  use,  the  conductivity  cell  was  filled  with  Triton  for  one  hour 
and  then  drained.  Triton  (Octyl  Phenol  Ethoxylate)  is  a  mild  non-ionic  surfactant  (Sea- 
Bird  Electronics  2008).  Flushing  of  the  conductivity  cell  prior  to  and  after  profiling  was 
accomplished  in  accordance  with  manufacturer’s  procedures.  When  the  CTD  was  not 
being  used,  the  inlet  to  the  conductivity  cell  and  the  outlet  of  the  SBE  5M  Mini 
Submersible  Pump  were  connected  with  a  length  of  Tygon  tubing  to  prevent  airborne 
contaminants  from  entering  the  conductivity  cell. 

When  not  in  use,  the  pH  probe  was  kept  clean  and  wet  in  its  buffer  solution 
saturated  with  potassium  chloride  (KC1).  This  soaker  fluid  is  acidic  (pH  4)  and  is 
designed  to  keep  the  electrode  free  of  marine  fouling.  The  internal  electrolyte  contains 
antibacterial  inhibitors  and  is  non-organic  to  support  use  in  the  ocean  (Sea-Bird 
Electronics  2010). 

The  pH  sensor  was  calibrated  daily  using  buffer  solutions  of  pH  10,  7,  and  4. 
Prior  to  calibration,  the  sensor  was  soaked  in  DI  water  for  30  minutes.  This  soaking 
allowed  the  charge  balance  around  the  sensor  to  stabilize  and  was  perfonned  in 
accordance  with  the  Sea-Bird  Electronics  calibration  procedures. 

E.  GPS  LOGGING 

The  position  of  the  profiles  was  logged  for  the  start  of  each  cast.  GPS  position 
was  logged  continuously  throughout  the  course  of  APLIS  2009.  For  this,  an  external 
antenna  was  mounted  to  the  top  of  the  working  hut  and  connected  to  a  computer  that 
logged  the  position  at  1Hz.  A  handheld  GPS  receiver  logged  the  positions  of  the  radial 
casts. 
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V.  ANALYTIC  METHOD  AND  DATA  ANALYSIS 


A.  DATA  PROCESSING 

The  MATLAB  functions  referenced  below  were  developed  by  Stanton’s  Ocean 
Turbulence  Group  (OTG)  at  the  Naval  Postgraduate  School  (Stanton  2003). 

1.  Profile  Isolation 

Usually,  there  were  many  profdes  together  in  one  data  fde  since  data  collection 
commenced  upon  activation  of  the  magnetic  reed  switch,  which  supplied  power  to  the  CTD. 
Many  profdes  were  conducted  before  removing  the  CTD  from  the  water  and  securing  power. 
After  the  CTD  data  was  converted  from  hexadecimal  to  ascii  format  using  the  Sea-Bird  Data 
Processing  program,  the  profdes  were  separated  out  for  comparison.  The  periods  during 
equilibration,  between  profdes,  and  before  securing  the  CTD  were  removed  from  the  profdes. 
These  data  were  saved  for  later  reference,  being  useful  for  studying  sensor  equilibration 
behavior  and  sensor  stability.  They  were,  however,  not  applicable  when  comparing  individual 
profdes. 

2.  Binning 

The  profde  data  were  consolidated  into  one-meter  bins  to  facilitate  processing  and 
comparison.  To  do  this,  the  OTG  MATLAB  function  ctd_bin.ni  was  used.  It  averaged  the 
values  between  z-0.5m  to  z+0.5m  to  be  the  value  for  that  profde  at  z  meters.  Profiling  at 
0.5m/s  gave  approximately  eight  measurements  to  average  for  each  binned  meter,  ctd_bin.ni 
also  separated  out  individual  profdes  (making  a  down-up  set  into  a  downward  profde  and  an 
upward  profde).  In  order  to  isolate  the  turning  points  between  a  down  and  up  profde,  a  third 
order  low-pass  Butterworth  filter  was  temporarily  applied  to  the  depth  vector.  This  allowed 
the  actual  profdes  to  be  separated  from  pressure  fluctuations  due  to  sway  or  uneven  motion 
due  to  operation  of  the  winch. 

3.  De-spiking  Measurements 

An  OTG  MATLAB  function  despike. m  was  used  to  de-spike  conductivity.  The 
most  probable  cause  of  spiking  is  biota  being  swept  into  the  suction  of  the  conductivity 
cell.  It  removed  spikes  exceeding  a  threshold  value  (40  mS/cm)  by  interpolating  between 
the  adjacent  points. 
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4.  Temperature  Structure 

The  temperature  in  the  mixed  layer,  the  upper  30m,  was  very  uniform.  In  the 
upper  ten  meters,  the  minimum  temperature  was  -1.45 1°C,  the  maximum  was  -1.402°C, 
with  a  mean  of  -1.430°C  and  standard  deviation  of  0.010  °C.  This  is  expected,  because  it 
is  very  close  to  the  freezing  temperature  of  the  seawater  at  the  salinity  measured  in  the 
mixed  layer.  The  salinities  ranged  from  a  minimum  of  25.753  to  a  maximum  of  26.751, 
with  a  mean  of  26.360  and  standard  deviation  of  0.103.  These  values  are  displayed  in 
Table  3.  Statistics  for  the  thermocline  and  halocline,  depth  of  greatest  temperature  and 
salinity  change,  respectively,  are  also  shown  in  Table  3.  A  plot  of  the  thennocline  and 
halocline  data  is  shown  in  Figure  10.  Ten  of  the  measurements  on  day  88  were  performed 
as  part  of  the  radial  surveys,  five  stops  at  9.4km  increments  from  APLIS.  Measured 
temperature  structure  is  shown  in  Figure  1 1 . 


min 

max 

F 

a 

temperature  (°C),  upper  10m 

-1.451 

-1.402 

-1.431 

0.006 

salinity,  upper  10m 

25.753 

26.751 

26.360 

0.103 

thermocline  depth  (m) 

29 

36 

32.190 

1.300 

halocline  depth  (m) 

29 

36 

32.174 

1.082 

Table  3.  Temperature  and  salinity  minimum,  maximum,  mean,  and  standard  deviation 
for  the  upper  10m.  Also  shown  are  the  statistics  for  the  thennocline  depth  and 

the  halocline  depth. 
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Thermocline  and  halocline,  APLIS  2009 


Figure  10.  Thermocline  and  halocline  depths  versus  time,  in  Ordinal  days.  The  green 
and  black  lines  represent  the  data  after  being  filtered  to  indicate  thermocline 

and  halocline  depth  trends. 

5.  Salinity  Calculation 

Salinity  was  calculated  using  the  MATLAB  function  salinity.m.  This  function 
uses  the  Practical  Salinity  Scale  of  1978  equations  from  Fofonoff  and  Millard,  1983  in 
the  UNESCO  papers.  This  program  uses  temperature  (°C,  ITS-90),  conductivity  (mS/cm), 
and  depth  (as  decibars)  to  generate  unitless  salinity  from  the  ratio  of  measured 
conductivity  to  reference  conductivity.  The  reference  is  the  conductivity  of  standard 
seawater  of  practical  salinity  35,  T=15°C,  and  atmospheric  pressure.  By  definition,  it  is 
equal  to  the  conductivity  of  a  reference  solution  of  potassium  chloride  (KC1)  at  the  same 
temperature  and  pressure  (Fofonoff  and  Millard  1983).  The  salinity  structure  measured  at 
APLIS  is  shown  in  Figure  1 1. 

The  structure  of  the  salinity  profile  shows  the  mixed  layer  extending  to  about  32m 
depth.  The  bulge  at  100m  is  the  center  of  the  upper  halocline  layer  (UHL),  and  is 
primarily  of  Pacific  Ocean  origin  (Bates  et  al.  2009).  The  slightly  concave-up  section 
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from  about  150m  to  230m  is  the  lower  halocline  layer  (LHL),  made  up  of  waters 
originating  from  the  Atlantic  Ocean  but  modified  during  transit  (Bates  et  al.  2009).  The 
salinity  profile  changes  very  little  in  the  rest  of  the  water  column.  Bates  describes  two 
additional  layers  of  Atlantic  origin,  the  Atlantic  Water  Layer,  and  the  Arctic  Ocean  Deep 
Water,  but  salinity  is  very  near  35  for  both. 


Temperature  APLIS  2009  Salinity  APLIS  2009 


Figure  11.  Temperature  and  salinity  measured  at  APLIS.  Radial  survey  data  is 

excluded. 

6.  pH  Structure 

In  order  to  measure  pH,  the  voltage  was  measured  across  the  glass  electrode  of  a 
Seabird  Electronics  SBE-18  pH  sensor.  This  output  voltage  and  the  measured 
temperature  were  used  by  Sea-Bird’s  SeaSoft  Data  Processing  Software  to  calculate  pH. 
Equation  (5.1)  shows  how  the  pH  value  is  calculated.  Vout  is  the  output  voltage  of  the 
sensor  (0-5  volts),  T  is  the  measured  temperature  (K),  slope  and  offset  are  calibration 
coefficients. 
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(5.1) 


tt  „  ,  Vout-  offset 

pH  —  7  +  ( - 7 - ) 

1.98416x10 ~4  *T* slope 

When  calibrated  at  the  manufacturer  in  February  2009,  the  values  were  slope  =  4.3350 

and  offset  =  2.5277  ,  when  the  temperature  was  20.8°C.  The  sensor  was  calibrated  daily  at 

temperatures  that  ranged  from  19.3°C  to  25.4°C. 

During  the  experiments,  the  pH  sensor  had  an  unexpected  temporal  drift 
following  power-on  and  subsequent  change  from  room  temperature  to  the  -1.4°C 
seawater,  as  shown  in  Figure  12.  The  indicated  pH  would  generally  rise  as  the  time  in  the 
water  increased.  When  the  CTD  was  removed  from  the  water  and  powered  off  (for  data 
upload  and  equipment  maintenance),  and  then  lowered  again  into  the  water,  the  measured 
pH  was  lower. 


APLIS  2009  pH  as  indicated 


Figure  12.  pH  as  indicated.  Note  that  the  shape  of  the  structure  is  very  similar  for  the 

profiles. 
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This  drift  occurred  despite  allowing  the  CTD  to  remain  in  the  water  for  an 
average  of  eight  minutes  during  equilibration.  Following  a  step  change  in  pH,  the  SBE  18 
is  designed  to  respond  to  63%  of  the  final  pH  value  within  one  second.  Subsequent 
discussions  with  technical  staff  at  Sea-Bird  Electronics  showed  that  a  similar  effect  had 
been  seen  in  the  lab.  Taking  the  sensor  from  room  temperature  to  a  1°C  test  bath  during 
thermal  shock  tests,  the  manufacturer  noted  that  it  often  takes  up  to  15  minutes  for  the  pH 
sensor  to  reach  99%  of  the  final  value.  The  pH  sensor  used  for  this  research,  however, 
seemed  to  be  equilibrating  over  the  course  of  hours.  This  behavior,  as  observed  in  the 
time  series  at  the  beginning  of  APLIS,  is  shown  in  Figure  13. 
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Figure  13.  Measured  pH  at  10  meters  vs.  time,  during  the  time  series  at  APLIS.  This 
plot  demonstrates  the  effect  of  equilibrating  over  hours  instead  of  seconds  to 

minutes. 


Each  day  the  pH  sensor  was  calibrated  to  a  set  of  buffer  solutions  of  pH  10,  pH  7, 
and  pH  4  at  room  temperature.  The  indicated  pH  was  consistently  higher  than  the  pH  7 
and  pH  10  buffers,  and  lower  than  the  pH  4  buffer.  With  respect  to  the  pH  7  buffer,  the 
sensor  indicated  high  each  time  by  about  ApH  =  0. 1 187 . 
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The  time-dependent  behavior  was  not  observed  during  the  calibrations;  instead, 
the  indicated  pH  rapidly  approached  the  buffer  value.  It  was  noted  that  the  temporal  drift 
could  have  been  a  result  of  the  large  difference  between  the  temperature  of  the  air  in  the 
hut  and  the  temperature  of  the  seawater.  The  matter  is  being  pursued  with  the 
manufacturer,  but  it  could  have  been  an  issue  with  ionic  balance  between  the  buffers 
(which  have  a  deionized  water  base,  and  a  correspondingly  low  ionic  strength)  and  the 
seawater  (which  has  a  high  ionic  strength)  (see  Zeebe,  R.  E.  and  D.  A.  Wolf-Gladrow 
2001).  It  could  have  also  been  an  issue  with  the  ionic  balance  between  the  pH  4 
potassium  chloride  soaker  solution  that  the  pH  sensor  stayed  in  when  not  in  use,  and  the 
seawater.  It  could  have  also  been  a  combination  of  the  two. 

To  help  resolve  the  slow  turn-on  drift  in  the  pH  sensor,  Leif  Anderson,  a  polar 
ocean  chemist  from  the  University  of  Gothenburg,  was  consulted.  It  was  determined  that 
it  would  be  safest  to  match  the  pH  data  at  the  depth  of  least  environmental  variability. 
This  was  the  depth  of  minimum  temperature,  181m.  This  approach  is  shown  in  Figure  14. 
Also,  in  this  figure  are  pH  measurements  obtained  by  Anderson  from  a  cruise  through  the 
same  region  of  the  Beaufort  Sea  in  Fall  2005.  These  measurements  were  conducted  by 
obtaining  water  samples  at  discrete  depths  and  conducting  chemical  analyses  to 
determine  the  pH. 
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Comparison  of  indicated  pH  from  APLIS  2009 
and  Beaufort  Sea  Cruise  from  L.  Anderson  2005 


pH 

Figure  14.  The  red  curve  is  the  mean  of  two  profiles  from  Leif  Anderson’s  Fall  2005 
Beaufort  Cruise.  The  light  green  is  the  mean  of  the  collapsed  APLIS  data. 

Since  the  majority  of  pH  profiles  were  conducted  to  233m,  the  mean  of  the  4 
deeper  casts  was  aligned  with  the  end  of  the  shallower  profiles  to  extend  each  profile  to 
518m.  Due  to  low  diffusivity  and  slow  vertical  transport,  CCA  variability  is  largely 
confined  to  the  upper  hundreds  of  meters  (Bates  et  al.  2006)  and  is  largely  invariant  at 
depth  (Anderson  2005).  Because  pH  values  for  the  full  water  column  were  needed  to 
facilitate  absorption  calculations  along  the  ray  path,  a  decision  was  made  to  utilize 
Anderson’s  pH  data.  This  was  accomplished  by  maintaining  the  measured  pH  structure 
from  APLIS  and  translating  the  profiles  so  that  the  mean  of  the  measured  pH  profiles 
aligned  with  the  mean  of  the  Anderson’s  data.  The  transition  points  were  then  smoothed 
using  a  rectangular  convolution  filter.  The  resultant  pH  profiles  are  shown  in  Figure  15. 
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APLIS  2009, 

pHsws=-log([H]+[HS04]+[HF]) 


pH 

Figure  15.  Full  water  column  pH  profiles.  Shallow  data  is  from  APLIS  2009,  aligned 
with  deep  data  from  Anderson  2005,  using  the  Seawater  pH  Scale. 

There  were  also  other  choices  regarding  the  pH  scale  and  its  temperature- 
correction.  One  convention  is  to  adjust  pH  measurements  based  on  temperature  such  that 
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a  neutral  solution  will  be  measured  pH  7  at  25  °C.  The  convention  used  by  Anderson  was 
to  temperature-correct  pH  to  15°C.  Regarding  pH  scale,  there  are  different  scales  used  in 
chemical  oceanography  depending  on  the  focus  of  the  study  (Bates  1975).  The  scale  used 
by  Anderson  was  the  seawater  scale,  where  the  pH  (denoted  pHSws)  includes  the 
concentration  of  hydrogen  sulfate  ions  and  fluoride  ions,  as  well  as  hydrogen  ions,  as 
shown  in  equation(5.2).  [//  ]F  is  the  concentration  of  “free”  hydrogen  ions  (which 

includes  hydrate  complexes  like  ll:f)  ,  where  a  hydrogen  ion  combined  with  a  water 
molecule  to  form  hydronium);  [  HS04  ]  i  s  the  concentration  of  sulfate  ions;  and  [Iff]  is 

the  concentration  of  fluoride  ions.  The  fluoride  effect  is  small  since  the  concentration  of 
fluoride  ions  is  much  lower  than  that  of  the  sulfate  ions  (Zeebe  and  Wolf-Gladrow  2001). 

pHsws=-log([H+]F  +  [HSO-4]  +  [HF])  =  -\og[H+]sws  (5.2) 

If  the  pH  data  were  left  aligned  with  the  deep  Arctic  Ocean  pH  from  Anderson  to 
deal  with  the  temporal  drift,  then  the  APLIS  data  would  end  up  being  temperature 
corrected  to  15  °C  and  using  the  seawater  pH  scale.  In  the  field  of  acoustics,  however,  the 
free  pH  scale  is  used  ( pH  =  pHF  ),  which  accounts  for  the  concentration  of  hydrogen 
ions  alone,  as  shown  in  equation  (5.3)  (Zeebe  and  Wolf-Gladrow  2001). 

pHF=-\og[H+]F  (5.3) 

Because  of  this  difference  in  scales,  it  was  necessary  to  attempt  to  convert  the  pH 
to  the  free  scale.  pH  is  a  complex  measurement  variable  that  depends  on  the  relative 
concentrations  of  numerous  constituents  of  an  acid-base  system.  There  is  not  a 
straightforward  conversion  that  can  be  applied.  There  is,  however,  a  relation  that  shows 
that  the  difference  between  ( pHF - pHsws )  is  about  0.12  units  at  T=25  °C,  and  S=35. 
Applying  this  difference,  A pH  =  0.12  was  added  to  the  pH  profile  shown  in  Figure  15. 
The  final  pH  profiles  used  for  absorption  calculations  are  shown  in  Figure  16. 
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APLIS  2009 
pHfree=-log[H] 


Figure  16.  Free  scale  pH  profiles  used  for  absorption  calculations. 

B.  DATA  FOR  DEEP  WATER 

In  order  to  use  the  ray  tracing  method  for  acoustic  propagation,  the  sound  speed 
profile  for  the  entire  water  column  must  be  known,  and  this  is  determined  from  the  full 
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water  column  temperature  and  salinity  profiles.  As  discussed  in  the  previous  section,  the 
majority  of  data  collected  at  APLIS  was  to  a  depth  of  233m  with  four  deeper  profiles  to 
518m.  In  order  to  extend  the  depth  of  profiles  at  the  camp  through  the  whole  water 
column,  the  mean  of  the  four  deep  profiles  was  concatenated  with  and  aligned  to  the 
lower  15m  of  the  shallow  profiles.  The  measured  and  extended  profiles  for  temperature 
are  shown  in  Figure  17.  Similarly,  the  corresponding  plots  for  salinity  are  shown  in 
Figure  18. 


Measured  T  Measured  T,  extended  with  deep  mean 


Temperature  (C)  Temperature  (C) 

Figure  17.  (L)  Measured  temperature  profiles  (four  to  5 18m,  198  to  233m  average 

depth)  and  (R)  measured  temperature  profile  extended  to  5 1 8m  with  the  mean 

of  the  deep  profiles. 
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Measured  S  Measured  S,  extended  with  deep  mean 


Figure  18.  (L)  Measured  salinity  and  (R)  measured  salinity  extended  with  the  mean 

of  the  deep  profiles.  These  plots  are  nearly  identical,  highlighting  the  low 
variability  of  salinity  below  the  upper  ocean.  The  transition  regions  where  the 
shallow  profiles  are  extended  are  visible  when  comparing  the  areas  between 

200m  and  250m. 

The  next  step  for  extending  the  depth  of  the  T  and  S  profiles  was  to  incorporate 
data  collected  at  about  the  same  time  from  an  Ice-Tethered  Profiler  (ITP)  that  was  located 
near  APLIS.  The  Ice-Tethered  Profiler  data  were  collected  and  made  available  by  the  Ice- 
Tethered  Profiler  Program  based  at  the  Woods  Hole  Oceanographic  Institution 
(http://www.whoi.edu/itp).  Two  profilers  were  active  near  APLIS.  ITP  21  was  located 
454km  bearing  028°,  and  ITP  1 1  was  located  496km  bearing  057°.  The  data  from  ITP  21 
was  more  similar  to  the  deep  casts  from  APLIS  suggesting  a  more  similar  water  mass. 
Each  day  the  ITP  profiled  to  760m  and  measured  T  and  S.  Because  of  this  similarity 
between  APLIS  and  ITP  data,  the  mean  of  the  daily  ITP  21  data  from  500m  to  760m  was 
added  to  the  bottom  of  the  extended  APLIS  data,  assuming  that  the  mean  was  a  good 
representation  of  T  and  S  in  this  depth  range  over  this  time.  The  ITP  data  is  provided  with 
2m  bin  spacing.  It  was  linearly  interpolated  to  obtain  lm  bin  widths. 
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For  data  in  the  rest  of  the  water  column,  a  data  set  was  selected  from  the  World 
Ocean  Atlas  (WOA)  of  2001.  This  data  set  was  provided  by  the  National  Oceanographic 
Data  Center  (National  Oceanographic  Data  Center  2007).  The  WOA  profiles  provide  data 
at  33  standard  depth  levels  and  were  linearly  interpolated  to  obtain  one-meter  bin  widths. 
The  WOA  database  has  data  for  half-degree  increments  of  latitude  and  longitude. 
Selecting  the  closest  location  in  the  database  to  APLIS  (33km  North),  full  water  column 
profiles  for  T  and  S  were  obtained.  The  profiles  were  an  annual  average,  and  the  values 
from  760m  to  3500m  were  added  below  the  ITP  21  data.  Comparison  of  the  WOA,  ITP 
1 1,  ITP  21,  and  the  mean  of  the  deep  profiles  from  APLIS  are  shown  in  Figure  19. 


Compare  APLIS,  ITP  11  &  21,  WOA  2001  Compare  APLIS,  ITP  11  &  21,  WOA  2001 


Temperature  (C)  Salinity 


Figure  19.  (L)  Temperature  profiles  from  WOA2001  near  APLIS,  ITP  11,  ITP  21, 

and  the  mean  of  the  deep  profiles  obtained  at  APLIS.  (R)  Salinity  profiles 

from  the  same  sources. 

Combining  these  data  sources  allowed  for  the  estimation  of  full  water  column 
profiles  for  temperature  and  salinity,  while  capturing  the  measured  variability  in  the 
upper  ocean.  It  was  assumed  that  for  the  times  studied  the  variability  of  T  and  S  below 
760m  was  negligible,  and  the  annual  average  at  this  location  was  a  good  representation  of 
the  deep  ocean  T  and  S.  The  full  water  column  profiles  used  for  ray  path  and  absorption 
calculations  are  shown  in  Figure  20. 
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Combined  Temperature  Profiles  Combined  Salinity  Profiles 


Figure  20.  Full  water  column  T  and  S  profiles  used  to  calculating  sound  speed  and 

absorption. 


C.  SOUND  SPEED  PROFILES 

1.  Sound  Speed  Calculation 

Using  the  inputs  of  temperature,  salinity,  and  depth,  from  the  full  water  column 
profiles  allowed  the  calculation  of  sound  speed  profiles  for  the  full  water  column,  as 
shown  in  Figure  20.  The  Woods  Hole  Oceanographic  Institute  Acoustic  Group’s 
MATLAB  Oceanography  Utilities  function  soundspeed.m  was  used  to  calculate  the 
sound  speed  profiles.  It  can  implement  several  sound  speed  equation  sets.  The  algorithm 
used  for  ray  tracing  and  absorption  calculations  was  CML  1 .0,  the  equations  by  Chen  and 
Millero  in  1977,  as  modified  by  Millero  and  Li  in  1994.  This  algorithm  is  the  U.S.  Navy 
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Standard  Sound  Speed  Algorithm  and  was  used  for  this  study.  It  was  selected  due  to  best 
performance  in  extensive  testing  and  comparison  with  other  algorithms.  The  current, 
effective  January  26,  2010,  Navy  Standard  is  CML  1.1,  which  is  the  same  equation  set, 
but  with  expanded  limits  for  the  T,  S,  and  D  inputs  (Moskal  2010).  A  comparison  of 
sound  speed  profiles  calculated  by  Mackenzie,  1960;  Leroy  et  al.  2008;  Del  Grasso, 
1974;  and  Chen-Millero-Li,  1994  is  shown  in  Figure  21.  The  plot  on  the  left  was 
generated  using  the  mean  profiles  for  temperature  and  salinity  from  APLIS  2009.  The 
plot  on  the  right  used  all  of  the  temperature  and  salinity  profiles  to  calculate  SSPs  for  all 
of  the  equation  sets.  It  is  important  to  note  that  while  some  data  was  obtained  to  518m, 
98%  of  the  profiles  were  to  a  mean  depth  of  233m.  This  figure  illustrates  that  there  is 
little  spread  in  the  shallow  water.  While  the  profiles  diverge  below  1000m,  there  is  less 
than  1.2m/s  (0.08%)  difference  at  3500m.  The  same  data,  zoomed  in  to  the  upper  250m, 
is  shown  in  Figure  22.  This  illustrates  that  given  the  T  and  S  profiles  obtained  at  APLIS 
2009,  the  spread  in  the  SSP  in  the  upper  ocean  due  to  the  selected  algorithm  is  dominated 
by  the  spread  in  the  SSP  due  to  variability  in  the  T  and  S  profiles. 
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SSPs  from  Mean  APLIS  2009  data  SSPs  for  all  APLIS  2009  data 


Figure  2 1 .  CML  sound  speed  profde  (SSP)  in  black  compared  to  SSPs  calculated  by 

three  other  equation  sets.  The  subplot  on  the  left  shows  the  SSPs  calculated 
from  the  mean  T  and  S.  The  subplot  on  the  right  shows  all  of  the  SSPs,  for  all 
of  the  equation  sets,  generated  from  all  the  T  and  S  profiles. 
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SSPs  from  Mean  APLIS  2009  data  SSPs  for  all  APLIS  2009  data 


Figure  22.  CML  sound  speed  profde  (SSP)  in  black  compared  to  SSPs  calculated  by 
three  other  equation  sets.  The  subplot  on  the  left  shows  the  SSPs  calculated 
from  the  mean  T  and  S.  The  subplot  on  the  right  shows  the  all  of  the  SSPs,  for 
all  of  the  equation  sets,  generated  from  all  the  T  and  S  profiles  zoomed  in  to 

the  upper  250m. 

2.  Sound  Speed  Profile  Features 

These  figures  highlight  the  strongly  upward-refracting  environment  common  to 
the  Arctic  Ocean  and  the  fact  that  the  environment  is  most  strongly  upward-refracting  in 
the  upper  400m.  They  also  show  the  100m  deep  sound  channel  centered  at  150m.  Using 
the  mean  SSP,  the  channel  extended  from  89m  to  195m.  The  value  of  sound  speed  at  the 
top  of  the  channel  is  1442.77m/s,  and  the  sound  speed  minimum  was  1441.61m/s  at 
123m.  This  channel  was  a  key  factor  for  ray  propagation,  particularly  for  cases  when  the 
sound  source  was  located  at  150m. 
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D.  RAY  MODEL 

The  ray  model  used  was  based  on  an  Eigenray  model  developed  by  Prof.  John 
Colosi  and  modified  by  Ben  Jones  from  the  Naval  Postgraduate  School.  It  assumes  a 
perfectly  reflecting  surface  and  uses  the  image  method  to  model  ray  surface  reflections. 
In  order  to  sample  the  possible  eigenrays,  sixty  launch  angles  were  selected  for  the  model 
runs,  from  -15°  to  +15°.  Three  source  depths  (50m,  100m,  and  150m)  were  evaluated  in 
order  to  study  the  effects  of  source  depth  location  with  respect  to  the  sound  speed  channel 
and  other  structure  in  the  upper  ocean.  For  a  given  launch  angle  and  source  depth,  the 
model  detennined  the  sound  speed  and  sound  speed  gradient  at  that  point.  It  then  used  the 
ordinary  differential  equation  solver  MATLAB  function  ode45.m  to  solve  the  eikonal 
equation,  calculating  the  sound  speed  gradient  and  ray  trajectory  to  detennine  the  ray 
path  in  5m  intervals  to  a  maximum  range  of  50km.  The  ray  path  consists  of  a  depth  for 
each  horizontal  range  step  of  5m,  and  then  the  Pythagorean  Theorem  determines  the 
actual  path  length  in  that  range  step. 

1.  Inputs 

a.  Temperature  and  Salinity 

The  ray  tracing  was  conducted  for  each  profile  of  T  and  S.  For  each  set  of 
T  and  S  profiles,  the  sound  speed  profile  was  detennined  using  soundspeed.m  as 
described  in  the  previous  section.  The  T  and  S  profiles  were  also  used  for  the  calculation 
of  absorption  as  will  be  discussed  in  the  next  section. 

b.  pH 

The  ray  tracing  is  not  dependent  on  pH.  pH  was  included  for  its 
contribution  to  absorption  and  transmission  loss,  which  was  implemented  as  described 
above.  pH  structure  is  generally  not  considered  in  absorption  calculations.  When 
absorption  is  calculated  with  constant  pH,  a  single  value  is  used  for  the  whole  water 
column. 

In  order  to  study  the  effects  of  pH  structure,  as  well  as  the  effects  of 
lowering  pH  due  to  CCfi  uptake,  a  pH  effect  switch  was  included  in  the  transmission  loss 
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portion  of  the  ray  model.  To  detennine  the  effects  of  the  pH  structure,  the  model  ran  once 
for  each  case  with  pH  assumed  to  be  constant,  and  once  with  the  measured  pH  profde  for 
the  whole  water  column. 

2.  Sound  Speed  Profile  Smoothing 

Since  the  profiles  were  processed  with  lm  bin  widths,  there  were  instances  where 
the  temperature  and  salinity  profiles  produced  sharp  corners  in  the  sound  speed  profile. 
When  these  comers  were  too  sharp,  the  differential  equation  solver  broke  down.  Whether 
a  profile  would  support  processing  with  ode45.m  was  dependent  on  the  source  depth  also. 
Fewer  profiles  ran  into  difficulties  with  a  source  at  150m.  In  order  to  have  most  of  the 
profiles  yield  solutions  while  still  preserving  the  structure  and  uniqueness  of  the  profiles, 
the  SSPs  were  low-pass  filtered  using  a  second  order  symmetric  Butterworth  digital  filter 
using  an  OTG  MATLAB  routine  Ipfilt.m.  This  served  to  filter  the  SSP  with  cutoff 
frequency  0.1  m1,  yielding  a  vertical  resolution  of  20m,  and  allowing  95%  of  the  cases  to 
yield  ray  traces. 

3.  Ray  Traces 

A  sample  sound  speed  profile  is  shown  in  Figure  23.  A  set  of  ray  traces  for  this 
SSP  is  shown  in  Figure  24.  For  this  figure,  the  source  is  at  150m,  near  the  axis  of  the 
sound  channel.  This  highlights  the  large  number  of  interactions  with  the  surface  (ice),  and 
shows  no  bottom  interactions.  The  eight  (of  60)  rays  with  the  narrowest  launch  angles 
(-1.78°  to  +1.78°)  remained  trapped  in  the  sound  channel  and  had  no  interactions  with  the 
surface.  For  the  same  SSP,  a  set  of  ray  traces  with  the  source  located  at  100m,  in  the 
sound  channel  but  near  the  sound  speed  maximum,  is  shown  in  Figure  25.  The  sound 
speed  maximum  for  this  profile  actually  occurs  at  87m.  In  this  case,  only  the  four 
narrowest  launch  angle  rays  (-0.76°  to  +0.76°)  were  trapped  in  the  sound  channel.  For  the 
same  conditions,  but  with  the  source  located  above  the  sound  channel  at  50m,  the  ray 
traces  are  shown  in  Figure  26.  This  figure  shows  the  much  higher  number  of  shallow- 
angle  surface  interactions.  In  this  case,  none  of  the  rays  was  trapped  in  the  sound  channel. 


48 


Figure  23. 


Sample  SSP,  Profile  3 
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Sound  speed  profde  that  was  used  for  the  following  ray  traces. 


Sample  Ray  Trace  (Profile  3.  f=1kHz,SD=150m)  from  APLIS  2009 
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Figure  24.  Set  of  ray  traces  for  APLIS  profde  #3,  source  depth  of  150m,  and 

frequency  of  1kHz.  Note  the  density  of  rays  in  the  sound  channel. 
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Sample  Ray  Trace  (Profile  3,  f=1KHz,  SD=100m)  from  APLIS  2009 


Range  (km) 

Figure  25.  Set  of  ray  traces  for  APLIS  profile  #3,  source  depth  of  100m,  and 
frequency  of  1kHz.  Note  that  there  are  fewer  rays  trapped  in  the  sound 
channel,  and  a  higher  density  of  rays  interacting  with  the  surface. 
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Sample  Ray  Trace  (Profile  3,  f=1KHz,  SD=50m)  from  APLIS  2009 


Figure  26.  Set  of  ray  traces  for  APLIS  profile  #3,  source  depth  of  50m,  and  frequency 
of  1kHz.  Note  that  there  are  no  rays  confined  to  the  sound  channel  and  the 
large  number  of  surface  interactions. 

The  rays  launched  from  50m,  above  the  sound  channel,  have  the  largest  number 
of  interactions  with  the  surface.  They  have  an  average  of  twice  as  many  surface 
reflections  compared  to  rays  launched  from  100m  or  150m.  This  will  factor  significantly 
into  the  structure  of  the  ice  scatter-absorption  contribution  to  transmission  loss,  which  is 
discussed  in  the  next  section.  The  number  of  surface  reflections  versus  launch  angle,  for 
each  source  depth,  is  shown  in  Figure  27.  and  shows  the  significance  of  the  position  with 
respect  to  the  sound  channel.  A  comparison  of  the  difference  in  structure  for  very  shallow 
angles  is  very  interesting.  The  rays  launched  from  50m  are  not  caught  in  the  sound 
channel  and  have  up  to  twice  as  many  surface  reflections,  whereas  some  rays  launched 
from  100m  never  exit  the  sound  channel.  An  even  greater  number  of  rays  launched  from 
150m  remain  trapped  in  the  channel  over  the  50km  range  of  study,  and  never  interact 
with  the  surface. 
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Number  of  Ice  Reflections  as  a  Function  of  Source  Depth  and  Frequency 


Launch  Angle  (Degrees) 


Figure  27.  Number  of  reflections  from  the  ice  surface,  as  a  function  of  frequency  for 
the  three  source  depths.  Note  how  the  negative  launch  angles  tend  to  have  one 
less  reflection  than  the  corresponding  positive  launch  angles. 

E.  TRANSMISSION  LOSS 

Using  ray  theory  allowed  for  a  comparison  of  the  relative  magnitudes  of 
transmission  loss  factors  by  keeping  the  effects  of  different  contributors  separate.  This  is 
an  approximation  for  convenience  because  the  focus  of  this  thesis  is  the  effects  of 
environmental  variability.  This  model  detennines  the  TL  for  a  given  path  and  is  not  a 
specific  source-to-receiver  transmission  loss  calculation  that  would  include  all  possible 
paths  and  the  consideration  of  the  coherent  vs.  incoherent  contributions  to  acoustic 
intensity. 

A  function  was  added  to  the  ray  model,  which  calculates  the  Francois-Garrison 
absorption  coefficient  at  each  step  along  the  ray  path  (with  the  current  profile-dependent 
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and  depth-dependent  values  of  T,  S,  and  pH),  then  multiplies  this  coefficient  by  each 
corresponding  step  path  length  to  give  a  TL  due  to  absorption  for  that  step.  All  the  step 
absorption  values  are  then  summed  to  give  a  value  for  the  absorption  contribution  to  TL 
for  that  ray.  Similarly,  the  geometrical  spreading  was  calculated  using  the  integral  path 
length  of  each  ray.  The  number  of  surface  reflections  was  determined  for  each  ray,  and 
the  ice  scatter-absorption  effect  was  calculated  using  a  dB-loss-per-ice-reflection  model 
from  Diachok,  1976. 


1.  Spherical  Spreading 

The  temperature  and  salinity  profiles  drive  the  sound  speed  profile.  The  gradient 
of  the  sound  speed  profile  detennines  the  curvature  of  the  acoustic  rays.  The  path  length 
differences  resulting  from  the  structure  of  the  sound  speed  profile  lead  to  an  angular 
dependence  in  the  transmission  loss.  This  path  length  difference  and  the  associated  TL 
are  shown  in  Figure  28.  There  is  also  a  dependence  on  the  source  depth  location.  The 
position  of  the  sound  source,  with  respect  to  the  sound  speed  profile  structure,  influences 
how  long  the  ray  path  will  be,  as  shown  in  Figure  28.  Note,  that  for  the  source  at  50m, 
there  is  a  relatively  stable  TL  at  narrow  angles.  This  is  because  rays  launched  at  very 
narrow  angles  follow  similar  paths.  For  comparison,  the  center  plot  shows  data  for  a 
source  depth  of  100m.  The  depression  in  the  center  results  from  the  rays  with  very 
narrow  launch  angles  (<0.8°)  that  are  trapped  in  the  sound  channel  and  have  shorter  path 
lengths,  never  undergoing  surface  reflections,  and  correspondingly  have  lower  TLgeo .  This 

phenomenon  is  seen  to  a  greater  degree  in  the  150m  source  depth  plot  on  the  right.  Since 
the  source  is  located  near  the  center  of  the  sound  channel,  rays  with  launch  angle  twice  as 
large  (<1.8°)  as  the  100m  case  are  trapped  in  the  sound  channel.  This  will  be  important  in 
the  next  section  when  comparing  the  10kHz  absorption  curves  at  150m.  These  angular 
confinements  agree  with  Snell’s  Law,  equation(5.4),  requirements,  c.  is  the  sound  speed 
and  0_  is  the  angle  of  the  ray  trajectory,  with  respect  to  the  horizontal,  at  depth  z  . 


Cl  _  C2 

cos  6X  cos  02 


(5.4) 
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From  Snell’s  Law  and  the  mean  sound  speed  profile  from  APLIS,  rays  launched 
from  100m  at  angles  <1.2°  will  be  trapped,  and  rays  launched  from  150m  at  angles  <2.2° 
will  be  trapped  in  the  sound  channel. 

Also  interesting  in  this  figure  is  the  positive/negative  launch  angle  effect.  Rays 
launched  with  a  positive  angle  tend  to  follow  a  path  similar  to  the  corresponding  negative 
angle,  but  have  an  additional  path  length  from  traveling  to  the  surface  and  back  to  the 
source  depth.  This  is  illustrated  by  the  points  on  the  left  side  of  the  curves  being  slightly 
lower  than  the  corresponding  points  on  the  right  of  the  curves  in  Figure  28.  Similarly,  the 
left-right  difference  due  to  the  positive  launch  angles  having  one  extra  surface  interaction 
was  also  seen  in  Figure  27. 


Source  depth  50m 


50 


94.1  50.8 


93.9  50  - l 


Source  depth  100m 

“i - 1 - 1 - 1 - 


Source  depth  150m 


93.9  50 


1  I  f\ 

-  -  - 

_ 

-15  -10  -5  0  5  10  15  -15  -10  -5  0  5  10  15  -15  -10  -5  0  5  10  15 

Launch  Angle  (Degrees)  Launch  Angle  (Degrees)  Launch  Angle  (Degrees) 


Figure  28.  For  one  profile,  plots  of  path  length  (blue)  and  spherical  TL  (green)  for 
source  depth  of  50m  (left),  100m  (center),  and  150m  (right). 

2.  Absorption  in  Seawater 

The  absorption  of  acoustic  energy  in  seawater  is  a  function  of  T,  S,  pH,  pressure, 
and  frequency.  The  absorption  profiles  for  the  three  frequencies  studied  are  shown  in 
Figure  29.  This  plot  was  created  using  the  Francois-Garrison  absorption  equation  (3.3), 
using  mean  T,  S,  and  pH  profiles.  The  mean  T  profile  was  obtained  by  averaging  the  T  at 
each  depth,  for  all  the  profiles  from  APLIS.  The  mean  profiles  for  S  and  pH  were 
calculated  in  the  same  way.  The  colored  curves  include  the  mean  profiles  as  discussed. 
The  black  curves  were  developed  using  the  mean  T  and  S  profiles,  but  with  a  single  pH 
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TLgeo  (dB) 


value  independent  of  depth.  This  value,  pH=7.98,  is  the  mean  of  all  the  pH  data,  for  all 
depths.  Since  pH  is  typically  not  considered  as  a  function  of  depth,  this  plot  illustrates  the 
differences  in  the  absorption  coefficient  when  pH  is  taken  to  be  a  single  value  (black 
curves),  and  when  the  measured  pH  is  used  (which  has  noticeable  depth-dependence; 
color  curves). 


Absorption  Profile,  including  measured  pH  vs.  constant  pH=7.98 


Figure  29.  Absorption  Profiles.  Color  curves  include  measured  parameters.  Black 

curves  consider  pH=7.98,  a  single  value  for  all  depths. 

This  figure  is  interesting  because  it  shows  how  the  relative  effects  of  the  T,  S,  pH 
and  frequency  influence  the  absorption  as  a  function  of  depth.  Locations  of  interest  on 
these  curves  are  summarized  in  Table  4. 


55 


Constant  pH 

Measured  pH 

100Hz 

lKHz 

lOKHz 

100Hz 

lKHz 

lOKHz 

Depth  (m) 

Absorption  Maximum 

1 

264 

231 

59 

64 

248 

Absorption  Minimum 
(upper  500m) 

455 

3 

3 

189 

185 

3 

Absorption  Minimum 

3500 

3500 

3500 

3500 

3500 

3500 

Table  4.  Absorption  Profile  points  of  interest.  Note  that  the  100Hz  and  1kHz  curves 
have  an  absorption  maximum  above  the  sound  channel  and  a  local  minimum 
in  the  channel,  while  the  10kHz  curve  has  a  maximum  absorption  coefficient 
below  the  channel  and  a  local  minimum  near  the  surface. 

The  effects  of  pH  structure  are  illustrated  in  the  next  six  figures.  The  surface  is 
assumed  to  be  perfectly  reflecting.  Note  that  geometrical  spreading  accounts  for  about 
94dB  in  these  plots.  The  upper  row  represents  model  runs  with  a  single  pH  value, 
independent  of  depth.  The  lower  row  includes  the  structure  in  the  measured  pH  profiles. 
The  columns  represent  source  frequencies:  the  one  on  the  left  for  100Hz,  the  center 
column  for  1kHz,  and  the  column  on  the  right  for  10kHz.  Note  the  scales  on  the  left  of 
each  plot  and  that  the  overall  magnitude  difference  is  small  when  considering  the  effects 
of  pH  structure.  The  thin  black  curve  is  the  mean  of  the  TL  data  in  each  plot,  and  the 
dashed  black  line  is  +  1  standard  deviation  from  the  mean.  This  description  applies  to 
Figure  30.  through  Figure  32.  The  transmission  loss  versus  launch  angle,  for  all  the 
profiles  at  APLIS  (i.e.  excluding  the  radial  survey),  are  shown  for  rays  launched  from 
50m  in  Figure  30.  Similar  plots  for  source  depth  of  100m  and  150m  are  shown  in  Figure 
31  and  32. 
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TL  (dB)  TL  (dB) 


f=0.1  kHz,  SourceDepth=50m,  pH  constant  f=1  kHz,  SourceDepth=50m,  pH  constant  f=  1 0  kHz,  SourceDepth=50m,  pH  constant 


Figure  30.  Transmission  loss  vs.  Launch  Angle.  Source  depth  =  50m. 
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TL  (dB)  TL  (dB) 


f=0.1  kHz,  SourceDepth=100m,  pH  constant  f=1  kHz,  SourceDepth=100m,  pH  constant  f=10  kHz,  SourceDepth=100m,  pH  constant 


f=0.1  kHz,  SourceDepth=100m,  pH  data  f^1  kHz,  SourceDepth=100m,  pH  data  f=10  kHz,  SourceDepth=100m,  pH  data 


Figure  3 1 .  Transmission  loss  vs.  Launch  Angle.  Source  depth  =  100m. 


58 


f=0.1  kHz,  SourceDepth=150m,  pH  constant  f=1  kHz,  SourceDepth=150m,  pH  constant  f=10  kHz,  SourceDepth=150m,  pH  constant 


Launch  Angle  (Degrees) 


Launch  Angle  (Degrees) 


Launch  Angle  (Degrees) 


f=0.1  kHz,  SourceDepth=150m,  pH  data  f=1  kHz,  SourceDepth=150m,  pH  data 


Launch  Angle  (Degrees)  Launch  Angle  (Degrees) 


f=10  kHz,  SourceDepth=150m,  pH  data 


Figure  32.  Transmission  loss  vs.  Launch  Angle.  Source  depth  =  150m. 


The  top-left  and  top-middle  plots  of  the  TL  vs.  Angle  plots  are  very  similar  to  the 
corresponding  plots  showing  the  contribution  of  spherical  spreading  in  Figure  28. 
Comparing  the  50m  source  depth  plots  of  Figure  30.  to  that  of  Figure  28.  the  100Hz  plot 
is  about  1.07  dB  higher,  and  the  1kHz  plot  is  about  3  dB  higher  as  a  result  of  absorption. 
The  10kHz  plot  on  the  top-right  of  Figure  30  is  significantly  different  from  the  TL  due  to 
spherical  spreading.  This  results  partially  from  the  increased  absorption  at  higher 
frequency.  At  medium  (5°-10°)  launch  angles;  there  is  higher  absorption  due  to  more  time 
spent  in  areas  with  a  larger  absorption  coefficient,  as  shown  in  Figure  29.  The  narrow 
launch  angles  give  more  time  close  to  the  surface,  where  there  is  the  absorption  minimum 

for  10kHz.  There  is  significant  variability  in  the  TL  for  narrow  angles  at  this  frequency 
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resulting  from  the  structure  of  the  absorption  profile  and  the  higher  magnitude  of  the 
absorption  coefficients  at  higher  frequency.  Also  at  10kHz,  the  steep  angles  have  the 
longest  ray  paths,  but  are  near  the  lowest  TL  because  they  also  spend  more  time  in  deeper 
waters  where  the  absorption  coefficients  are  again  low.  There  is  a  similar  structure 
noticeable  in  the  1kHz  plot  (in  the  top  row,  considering  a  single  value  for  pH),  but  with  a 
significantly  lower  amplitude,  resulting  from  the  lower  TL  at  lower  frequency.  This  effect 
would  be  opposite  for  100Hz,  since  the  pH  constant  absorption  profile  has  a  maximum  at 
the  surface,  but  this  is  negligible  because  the  absorption  is  so  small  at  low  frequency. 

These  absorption  profile  effects  are  accentuated  for  Figures  31  and  32,  when  the 
sound  source  is  located  within  the  sound  channel  where  the  pH  profile  has  more 
variability.  This  is  illustrated  by  comparing  the  50m  source  depth  in  Figure  30  with  the 
100m  source  depth  in  Figure  31.  Again,  the  wider  range  of  launch  angles  being  trapped  in 
the  sound  channel  is  shown  for  the  150m  source  depth  plots  of  Figure  32.  The  absorption 
profile  for  10kHz  has  a  larger  contribution  to  TL  when  the  rays  are  near  the  sound 
channel,  located  near  the  maximum  absorption  coefficient  at  248m  in  Figure  29.  There  is 
the  opposite  effect  for  the  100Hz  and  1kHz  plots,  since  they  have  minimum  absorption 
coefficients  near  the  sound  channel  axis,  and  the  reduced  TL  for  narrow  angles  is  evident 
in  the  bottom-left  and  bottom-center  plots  in  Figures  3 1  and  32. 

In  order  to  see  more  clearly  the  increased  spread  in  the  absorption  component  of 
TL  as  a  result  of  pH  profile  structure,  histograms  were  created  which  count  the  number  of 
occurrences  of  a  TL  of  a  given  magnitude.  These  are  shown  in  Figures  33-35,  for  each 
source  depth,  frequency  and  pH  case.  These  figures  also  illustrate  the  increased  spread  of 
absorption  that  the  measured  pH  structure  has  on  the  different  frequencies.  The  relative 
spread  is  much  larger  at  1kHz  than  the  other  frequencies,  and  that  the  effect  at  100Hz  is 
greater  than  that  at  10kHz.  At  10kHz,  the  most  noticeable  change  is  that  the  peak  near 
155  dB  is  lowered  by  about  1  dB  when  the  measured  pH  profile  is  included,  but  the 
spread  is  similar. 
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#  of  Events  #  of  Events 


1=0. 1  kHz,  SourceDepth=50m,  pH  constant  1=1  kHz,  SourceDepth=50m,  pH  constant  f=10  kHz,  SourceDepth=50m,  pH  constant 
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Figure  33.  Histograms  illustrating  the  spread  of  TL,  for  a  perfectly  reflecting  surface, 
with  Source  Depth  =50m.  These  correspond  to  the  plots  of  Figure  30. 
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#  of  Events  #  of  Events 


f=0.1  kHz,  SourceDepth=100m,  pH  constant  f=1  kHz,  SourceDepth=100m,  pH  constant 


f=10  kHz,  SourceDepth=100m,  pH  constant 


f=0.1  kHz,  SourceDepth=100m,  pH  measured  f=1  kHz,  SourceDepth=100m,  pH  measured 


f=10  kHz,  SourceDepth=100m,  pH  measured 


Figure  34.  Histograms  illustrating  the  spread  of  TL,  for  a  perfectly  reflecting  surface, 
with  Source  Depth  =100m.  These  correspond  to  the  plots  of  Figure  3 1 . 


f=0.1  kHz,  SourceDepth=150m,  pH  constant 


f=1  kHz,  SourceDepth=150m,  pH  constant  f=  1 0  kHz,  SourceDepth=150m,  pH  constant 


f=0.1  kHz,  SourceDepth=150m,  pH  measured  f=1  kHz,  SourceDepth=150m,  pH  measured  f=  1 0  kHz,  SourceDepth=150m,  pH  measured 


TL  (dB)  TL  (dB)  TL  (dB) 

Figure  35.  Histograms  illustrating  the  spread  of  TL,  for  a  perfectly  reflecting  surface, 
with  Source  Depth  =150m.  These  correspond  to  the  plots  of  Figure  32. 
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The  difference  is  small  between  the  magnitude  of  absorption  with  measured  pH 
and  the  absorption  with  a  single  pH  value.  While  some  spread  between  a  profile  with 
structure  and  a  profile  without  structure  was  noticeable  in  the  100Hz  and  10kHz  cases, 
the  largest  relative  change  was  seen  in  the  1kHz  cases.  The  magnitude  of  the  change  is 
from  about  0.1  dB  spread  for  the  single  pH  case  to  0.7-0. 8  dB  for  the  measured  pH  case. 
One  contribution  to  the  spread,  particularly  at  narrow  angles  is  a  large  pH  “channel”  that 
is  centered  at  181m,  and  is  about  160m  deep.  This  leads  to  reduced  absorption  for  rays 
traveling  through  this  channel  but  is  most  apparent  for  narrow  launch  angle  rays  from 
source  depths  in  the  pH  channel  (100m  and  150m).  It  is  important  to  note  that  these  cases 
still  assume  a  perfectly  reflecting  surface.  Scatter-absorption  effects  of  the  ice  are 
investigated  in  the  next  section. 

The  largest  relative  change  due  to  the  pH  structure  occurs  in  the  1kHz  plots, 
which  is  expected  due  to  the  relaxation  frequency  of  boric  acid  being  near  1kHz 
(dependent  on  T  and  S),  and  the  absorption  is  sensitive  to  changes  in  T,  S,  and  pH  (from 
equation  (3.6)).  Note  that  in  all  these  plots,  spherical  spreading  accounts  for  about  94dB 
of  transmission  loss.  Note  also  the  scale  for  each  of  the  plots.  The  scale  for  the  100Hz 
plots  covers  0.2dB,  while  the  1kHz  plots  cover  2.0dB,  and  the  10kHz  plots  cover  12dB. 
This  is  a  result  of  the  exponential  frequency  dependence  of  absorption  in  seawater  seen  in 
equation(3.3). 


3.  Ice  Scatter-Absorption  Effects 

a.  SCICEX  Data  Set  and  Diachok  Loss  Model 

Using  data  on  keel  depths  and  spatial  distributions  from  a  Science  Ice 
Exercise  (SCICEX)  submarine  cruise  conducted  in  the  Western  Arctic  Ocean  in  1988 
(NSIDC  1998),  and  equation  (3.7)  for  dB  lost  per  bounce  from  the  ice  (Diachok  1976), 
transmission  loss  due  to  ice  scatter-absorption  was  calculated.  It  contained  407 
measurements  of  keel  draft  and  the  distance  to  the  next  keel.  This  distance  was  inverted 
to  develop  a  spatial  frequency.  The  SCICEX  data  set  defines  a  keel  as  having  a  draft 
equal  to  or  exceeding  5m.  As  a  result  of  the  thinning  of  Arctic  sea  ice  and  the  defining  of 
a  keel  as  drafting  greater  than  5m,  this  dataset  from  1988  was  expected  to  yield  ice  drafts 
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that  are  deeper  than  average  for  2009.  Since  thinner  ice  can  ridge  more  easily,  the  spatial 
distribution  could  be  lower  than  average.  It  was  assumed  that  these  combined  effects 
would  mitigate  each  other  and  lead  to  order  of  magnitude  accuracy.  As  stated  in  Chapter 
III,  the  Diachok  high-frequency  model  used  is  independent  of  the  specific  frequency, 
provided  that  the  product  of  the  wavenumber  and  the  mean  keel  draft  is  greater  than 
unity,  kd>l,  which  is  the  case  for  all  frequencies  when  using  the  SCICEX  dataset. 

Each  time  a  ray  was  to  interact  with  the  ice,  a  matched  set  of  keel  draft  (d) 
and  spatial  frequency  (N)  was  selected  at  random  from  the  measured  distribution,  the 
incident  angle  was  determined,  and  the  reflectivity  was  calculated  following 
equation(3.7).  It  was  important  to  calculate  N  and  d  from  the  same  random  draw  since  the 
spatial  frequency  of  ridges  is  related  to  the  draft,  specifically  that  larger  ridges  tend  to 
have  larger  distances  between  them,  as  shown  in  Figure  36.  As  an  example,  reflectivity  as 
a  function  of  incident  angle  is  plotted  for  keel  spatial  frequency  N=9.5km‘l  and  keel  draft 
d=4.3m  in  Figure  37.  Also  shown  in  this  figure  is  the  SCICEX  reflectivity  curve 
calculated  from  the  mean  keel  draft  (d=8.5m),  and  the  mean  spatial  frequency  of  keels 
with  draft  between  8m  and  9m.  This  resulted  in  N=1 1.8km'1.  Figure  37  shows  that  the 
data  set  from  SCICEX  results  in  greater  loss  per  ice  reflection  than  Diachok’s  example 
parameters.  It  is  also  evident  from  equation(3.7)  and  this  figure  that  large  reflection 
losses  are  generated  near  6°. 
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Figure  36.  Spatial  Frequency  vs.  Keel  Draft  for  SCICEX  1988  data.  This  figure 
shows  that  the  spatial  frequency  is  higher  for  smaller  keels,  and  lower  for 

larger  keels. 
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Figure  37.  Comparison  of  Reflectivity  vs.  Angle  for  Diachok’s  N=9.5  km'1  and 

d=4.3m,  and  the  mean  ice  draft  data  from  SCICEX,  N=1 1.8  km'1  and  d=8.5m. 
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Incident  angles  that  were  observed  in  this  study  from  the  source  depths 
and  sound  speed  profdes  are  shown  in  Figure  38.  The  majority  of  incident  angles  have 
very  low  reflectivity  and  correspond  to  high  TL  due  to  ice  scatter-absorption.  These  low 
incident  angles  result  from  narrow  launch  angles.  The  observations  in  Figure  38  agree 
well  with  Snell’s  Law,  which,  for  the  mean  SSP  from  APLIS  dictates  a  possible  range  of 
incident  angles  from  6°-16°  for  50m  source  depth,  and  about  7°-17°  for  100m  and  150m 
source  depths. 
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Figure  38.  Observations  of  angles  for  rays  when  incident  upon  the  ice.  For  the  top 
plot,  Source  Depth  =  50m.  For  the  center  plot,  Source  Depth  =  100m.  For  the 
bottom  plot,  Source  Depth  =  150m.  Note  that  the  scale  on  the  left  is  different 
for  the  50m  case  because  of  the  much  higher  number  of  surface  reflections. 


b.  TL  Results  with  Ice  Scatter-Absorption  Included 

The  following  six  TL  vs.  Launch  Angle  figures  have  the  same  format  as 
the  previous  six.  Each  figure  represents  rays  launched  from  one  source  depth.  The  left- 
hand  column  is  for  a  frequency  of  100Hz,  the  center  column  is  for  a  frequency  of  1kHz, 
and  the  right-hand  column  is  for  a  frequency  of  lOkFIz.  The  top  row  considers  pFI  to  be 
constant,  and  the  bottom  row  uses  the  measured  pH  profiles  when  calculating  absorption. 
The  heavy  black  curve  is  the  mean  of  the  TL  when  ice  effects  are  considered,  and  the 
heavy  dashed  black  curves  are  +  1  standard  deviation  from  the  mean.  The  red  curves 
show  the  total  TL,  being  the  sum  contributions  of  spherical  spreading,  absorption,  and  ice 
scatter-absorption.  The  blue  curves  contain  the  same  data  as  the  corresponding  figures 
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when  the  surface  was  assumed  to  be  perfectly  reflecting  (Figures  30-32).  The  effects  of 
the  interactions  with  the  ice  now  dwarf  the  blue  curves. 


These  TL  vs.  Angle  plots  and  the  histograms  that  follow  are  also  useful  to 
see  that  the  ice  effect  does  not  follow  a  Gaussian  distribution.  This  was  as  expected, 
based  on  the  way  the  SCICEX  data  set  defined  a  keel  as  a  draft  exceeding  5m. 


f=0.1  kHz,  SourceDepth=50m,  pH  constant  f=1  kHz,  SourceDepth=50m,  pH  constant  f=10  kHz,  SourceDepth=50m,  pH  constant 


f=0. 1  kHz,  SourceDepth=50m,  pH  data 


f=1  kHz,  SourceDepth=50m,  pH  data 


f=10  kHz,  SourceDepth=50m,  pH  data 


Figure  39.  Transmission  Loss  vs.  Launch  Angle  with  ice  scatter-absorption  effects 
included.  Source  Depth  =  50m.  The  blue  curves  which  are  barely  visible, 
contain  the  same  data  as  in  Ligure  30.  Note  the  scale  on  the  vertical  axis, 
which  shows  the  huge  increase  in  TL  in  dB  due  to  ice  effects. 
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Figure  40.  Transmission  Loss  vs.  Launch  Angle  with  ice  scatter-absorption  effects 
included.  Source  Depth  =  100m.  Note  the  scale  on  the  vertical  axis,  compared 
to  the  previous  figure,  and  the  effects  of  rays  getting  confined  to  the  sound 

channel. 


69 


Launch  Angle  (Degrees) 
f=0. 1  kHz,  SourceDepth=150m,  pH  data 

500  I - 1 - i - i - 1 - 1 - 

iii 

iii 

400  -  J-  -  -  -i - 1 - 1 - 1 - 


-15  -10  -5  0  5  10  15 

Launch  Angle  (Degrees) 


Launch  Angle  (Degrees) 
f=1  kHz,  SourceDepth=150m,  pH  data 


500 

400 


-15  -10  -5  0  5  10  15 

Launch  Angle  (Degrees) 


Launch  Angle  (Degrees) 


-15  -10  -5  0  5  10  15 

Launch  Angle  (Degrees) 


Figure  4 1 .  Transmission  Loss  vs.  Launch  Angle  with  ice  scatter-absorption  effects 

included.  Source  Depth  =  150m. 

As  in  subsection  two.  above,  where  histograms  were  used  to  illustrate  the 
distribution  of  TL  for  the  different  cases,  histograms  of  the  TL  data  in  Figures  39^4 1  are 
shown  in  Figures  42^44.  These  histograms  show  again  that  the  predicted  spread  of  TL  for 
paths  that  involve  ice  interactions  dominates  what  little  variability  was  present  due  to  pFI 
considerations. 
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Figure  42.  Histograms  illustrating  the  spread  of  TL  when  ice  scatter-absorption  is 
included.  Source  Depth  =  50m.  These  histograms  correspond  to  the  plots  in 

Figure  39. 
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Figure  43.  Histograms  illustrating  the  spread  of  TL  when  ice  scatter-absorption  is 
included.  Source  Depth  =  100m.  These  histograms  correspond  to  the  plots  in 

Figure  40. 
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Figure  44.  Histograms  illustrating  the  spread  of  TL  when  ice  scatter-absorption  is 
included.  Source  Depth  =  150m.  These  histograms  correspond  to  the  plots  in 

Figure  41. 


4.  Relative  Magnitude  of  Effects 

When  the  surface  was  considered  to  be  perfectly  reflecting,  the  largest 
contribution  to  the  magnitude  of  transmission  loss,  over  the  50km  studied  and  the 
frequencies  considered,  was  the  spherical  geometrical  spreading.  It  accounted  for 
94dB/50km.  Being  only  dependent  on  range;  the  only  variability  was  due  to  very  small 
differences  in  ray  path  length.  Over  50km,  the  maximum  increase  in  path-length  is  about 
500m,  and  occurs  with  the  steepest  launch  angles.  This  results  in  an  increase  in  TL  of 
only  0.09dB,  one  of  the  smallest  sources  of  variability.  Of  course,  it  is  possible  that  a 
more  sophisticated  model  for  geometrical  transmission  loss  would  show  higher 
variability.  In  particular,  the  variability  due  to  uncertainty  in  the  exact  SSP  would  be 
interesting  to  quantify.  For  1kHz  and  10kHz,  the  next  largest  contribution  to  TL  is  the 
absorption  in  the  water  column.  At  1kHz,  this  adds  about  3dB  to  TL  when  a  single  depth- 
independent  pH  value  is  used,  and  adds  2-3dB  when  the  measured  pH  structure  is 
included.  The  variability  in  TL  due  to  absorption  at  1kHz  is  low,  however,  at  +0.5dB.  At 
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10kHz,  absorption  contributes  50-62dB  of  TL,  depending  on  launch  angle/ray  path.  The 
pH  effect  at  this  frequency  is  up  to  1  dB,  but  this  is  dominated  by  the  variability  due  to 
launch  angle.  At  100Hz,  the  absorption  contribution  to  TL  is  very  small  (O.ldB),  and  the 
variability  is  also  very  small  (<0.1dB). 

The  TL  due  to  interactions  with  the  ice,  using  the  Diachok  model,  is  not 
frequency  dependent.  Therefore,  the  magnitude  of  the  predicted  TL  contribution  due  to 
ice  effects  is  the  same  for  all  frequencies.  The  largest  variations  are  seen  when  comparing 
the  TL  due  to  ice  for  different  source  depths.  For  example,  Figure  39  shows  extremely 
high  TL  values  due  to  high  numbers  of  surface  interactions.  In  Figures  40  and  41,  the 
rays  from  narrow  launch  angles  have  reduced  ice  interactions  due  to  the  effects  of  the 
sound  channel,  or  having  no  ice  interactions  as  a  result  of  being  completely  trapped  in  the 
sound  channel. 

Developed  from  a  dB/bounce  model  that  factors  in  the  spatial  frequency  of  ice 
keels  and  the  distribution  of  keel  drafts,  the  variability  in  under-ice  structure  is  the  largest 
contributor  to  the  variability  in  transmission  loss  for  all  frequency  and  source  depth  cases. 
The  loss  due  to  the  ice  is  also  affected  by  the  incident  angle.  In  this  upward  refracting 
environment,  small  incident  angles  result  from  shallow  source  depths  and/or  narrow 
launch  angles.  The  combination  of  keel  distribution  and  incident  angles  plays  heavily  into 
the  amount  of  loss  due  to  the  ice  scatter-absorption,  as  shown  in  Figure  37. 

An  important  point  to  remember  is  that  despite  the  model  output  of  hundreds  of 
dB  of  TL  due  to  ice  effects,  a  ray  that  accumulates  a  very  high  TL  is  essentially 
attenuated  beyond  detection  at  that  range.  For  a  source  of  any  frequency  at  50m,  TL  due 
to  ice  interactions  increases  rapidly  corresponding  to  many  surface  bounces.  If  a  source  is 
shallow  and  close  to  the  ice,  steep  angles  and  relatively  smooth  ice  are  needed  for 
meaningful  amounts  of  energy  to  propagate. 

For  a  10kHz  signal,  the  cumulative  TL  is  very  high  regardless  of  source  depth. 
Only  the  rays  that  are  trapped  in  the  sound  channel  have  a  realistic  chance  to  contribute  at 
50km.  The  window  is  very  narrow  (±1.2°)  for  a  source  at  100m,  but  for  a  source  at 
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150m,  the  range  of  angles  that  can  be  trapped  is  greater  (+2.2°)  and  presumably  the 
corresponding  contribution  at  a  range  of  50km  would  be  greater. 

As  the  frequency  of  a  signal  goes  down,  the  overall  TL  goes  down  as  a  function 
of  the  frequency-dependent  absorption.  For  1kHz  and  100Hz  sources,  steep  launch  angles 
reduce  the  effects  of  ice  scatter-absorption,  but  the  lowest  TL  results  from  rays  confined 
to  the  sound  channel.  In  this  narrow  region  where  there  is  no  loss  due  to  ice  interactions, 
the  variability  due  to  consideration  of  pH  can  be  seen  the  center  and  left  plots  of  Figure 
31  and  32.  In  this  case,  a  difference  of  up  to  1  dB  can  be  seen  based  on  pH  vertical 
structure  and  profile  variability. 

F.  RADIAL  SURVEY 

It  had  been  assumed  that  the  sound  speed  profile  (SSP)  changes  only  slowly  as  a 
function  of  radial  distance,  and  that  the  SSPs  could  be  considered  range-independent.  The 
radial  survey  studied  this  hypothesis  to  test  the  validity  of  the  assumption. 

1.  Data 

For  the  radial  survey,  two  profiles  (one  down,  one  up)  were  obtained  at  each  of  5 
locations  over  50km  to  the  northeast.  Two  additional  profiles  were  obtained  upon  return 
to  the  camp  (the  0th  radial)  to  provide  a  maximum  spatial  separation  with  minimum 
temporal  separation.  Due  to  errors  with  the  CTD,  most  probably  freeze-up  of  the 
conductivity  cell,  three  down-profiles  were  discarded.  The  temperature  and  salinity  for 
these  profiles  were  erroneous  until  halfway  to  230m,  where  they  correlated  well  with  the 
corresponding  up-profiles.  Including  the  0th  radial,  there  were  nine  good  profiles  with 
which  to  investigate  radial  dependence  on  sound  speed  profile.  The  SSPs  generated  from 
these  profiles  are  shown  in  Figure  45. 

2.  Range  Dependence 

While  the  general  structure  of  the  sound  speed  profiles  is  as  expected,  there  is 
more  variability  than  anticipated.  Also  shown  in  Figure  45  is  the  mean  and  standard 
deviation  from  the  time  series,  which  consisted  of  124  profiles  over  36  hours.  During  this 

time,  the  camp  drifted  south  122m.  This  figure  illustrates  that  there  was  more  variability 
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over  7.3  hours  and  47km,  than  over  36  hours  and  0.1km,  suggesting  that  spatial 
variability  at  this  scale  dominates,  and  that  the  validity  of  the  range-independence 
assumption  should  be  investigated  further. 


Comparison  of  SSPs  for  6  profile  locations  with  mean  9.4km  separation 


Sound  Speed  (m/s) 


Figure  45.  Sound  speed  profdes  for  the  radial  survey.  The  light  green  curve  shows 
the  mean  of  the  SSPs  from  124  profiles  from  the  36-hour  time  series.  The 
black  dashed  curves  show  + 1  standard  deviation  from  the  mean.  This  figure 
shows  that  the  variability  in  sound  speed  profile  in  a  few  hours  over  47km  is 
greater  than  the  variability  in  the  SSP  over  36  hours  and  drifting  0. 1km. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY  OF  RESULTS 

T,  S,  pH,  and  sound  speed  profiles  were  input  to  a  ray  model  that  was  run  for 
several  scenarios,  varying  source  depth  and  frequency.  The  effects  of  pH  variability  were 
investigated  by  running  the  model  with  a  constant,  depth  independent  pH  for  the  whole 
water  column  and  comparing  that  to  model  runs  with  the  measured  pH  in  the  upper 
ocean.  Another  factor  of  variability  considered  was  the  effects  of  ice  interactions.  The 
relative  smoothness  or  roughness  of  the  ice  depends  on  the  frequency  of  sound  in 
question.  In  the  Diachok  model,  the  reflectivity  is  considered  to  be  independent  of 
frequency,  provided  that  the  frequency  and  the  mean  keel  draft  are  large  enough  (kd>  1 ). 
Diachok’ s  model  detennined  reflection  losses  as  a  function  of  incident  angle,  and  mean 
keel  draft  and  mean  spatial  distribution.  For  this  thesis,  a  dataset  of  keel  drafts  and 
spacing  was  randomly  sampled  in  order  to  provide  increased  accuracy  of  the  Diachok 
model.  These  effects  are  accounted  for  with  an  equation  giving  a  dB  lost  per  bounce  from 
the  ice.  The  ray  trace/absorption  model  was  used  to  look  at  the  fine-scale  variability  in 
the  measured  parameters  and  the  effects  on  transmission  loss  for  individual  rays  over  a 
range  of  50km.  Due  to  model  limitations,  the  geometrical  transmission  loss  for  each  path 
was  approximated  as  spherical  spreading  over  the  path  length. 

1.  No  Ice  Effects  /  Perfectly  Reflecting  Surface 

For  a  source  of  100Hz,  there  was  less  than  O.ldB  variability  (0.05dB  below  the 
mean  to  0.05dB  above  the  mean)  in  TL  over  50km,  for  all  calculated  sound  speed 
profiles,  regardless  of  pH  effects  or  source  depth  of  50m,  100m,  or  150m. 

For  a  source  at  1kHz  ,  when  pH  is  held  constant,  there  is  less  than  O.ldB 
variability  in  TL  regardless  of  source  depth.  When  the  variable  pH  from  the  data  is 
considered,  there  is  approximately  ldB  of  variability.  Sound  that  leaves  the  source  at 
steep  angles  has  less  interaction  with  this  variable  pH  field  and  consequently  lower 
variability.  At  shallower  angles,  and  for  a  source  depth  above  the  sound  channel  (50m), 
there  is  much  more  variability. 
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For  a  source  at  10kHz,  there  is  up  to  5dB  of  variability  in  TL  resulting  from  the 
structure  of  the  T,  S,  and  sound  speed  profdes.  The  difference  between  considering  pFI 
effects  or  not  yields  less  than  2dB  of  variability,  and  is  seen  primarily  in  narrow  angles, 
as  was  the  case  for  the  1kHz  source. 

2.  Ice  Scatter-Absorption  Included 

Considering  acoustic  energy  along  the  ray  path  lost  in  each  interaction  with  the 
ice,  using  Diachok’s  1976  equation,  there  was  the  most  significant  effect  when  the  source 
was  at  50m.  By  a  range  of  50km,  essentially  all  of  the  rays  were  attenuated  beyond 
detection.  This  was  expected  in  the  strongly  upward-refracting  environment,  and 
especially  since  the  sound  source  above  the  sound  channel  has  many  more  surface 
interactions  than  for  the  cases  when  the  source  is  in  the  channel.  Similarly,  the  effect  of 
including  ice  scatter-absorption  is  strongest  for  small  angles  that  generally  have  more 
surface  interactions.  The  only  rays  expected  to  contribute  significantly  at  50km  range  are 
those  from  a  narrow  angle,  for  a  source  of  frequency  less  than  about  lkFIz,  located  in  the 
sound  channel.  As  the  frequency  goes  down  below  about  70Hz  for  the  SCICEX  keel  draft 
data  (kd<l),  Diachok’s  high-frequency  equation  no  longer  applies.  Similarly,  for  a  source 
at  100Hz,  the  mean  keel  draft  would  have  to  be  shallower  than  2.3m  in  order  for 
Diachok’s  low-frequency  equation  to  apply.  This  relation  could  become  applicable 
because  as  the  ice  gets  thinner,  and  the  mean  keel  draft  depth  gets  smaller.  For  a  given 
frequency,  smoother  ice  reflects  sound  with  less  loss,  and  so  thinner  ice  and  shallower  ice 
ridges  can  lead  to  reduced  contribution  of  ice  scatter-absorption  at  increased  frequencies. 

B.  RECOMMENDATIONS 

1.  pH  and  SSP 

Arctic  maritime  activity  will  increase  corresponding  to  the  decline  in  sea  ice 
coverage.  Increased  CO2  in  the  atmosphere,  combined  with  reduced  ice  cover,  will  lead 
to  an  increased  uptake  of  CO?  in  the  upper  Arctic  Ocean.  This  will  contribute  to  a 
lowering  of  pH,  primarily  in  the  upper  ocean.  Absorption  of  acoustic  energy  will  increase 
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with  this  change  in  pH,  but  the  effect  is  small.  Uncertainty  in  transmission  loss  due  to 
sea  ice  scatter-absorption  is  expected  to  dominate  the  uncertainty  due  to  the  changing  pH 
structure. 

As  a  result  of  this  study,  continuous  pH  profding  is  deemed  unnecessary  for 
acoustical  applications  given  the  small  magnitude  of  pH  fluctuations  embedded  in  large 
sensor  drifts  and  the  operational  complexity  of  maintaining  accurate  calibration  of  the 
sensors.  Specifically  the  operational  monitoring  of  pH  by  submarines  is  discouraged. 
Above  1kHz,  the  variability  in  TL  due  to  SSP  variability  is  expected  to  exceed  that  of  pH 
effects,  especially  in  the  near-surface  environment.  Continued  study  of  temperature  and 
salinity  variability  in  the  upper  Arctic  Ocean  is  recommended. 

As  seen  in  Table  1.  and  Figures  30-35  and  39-44,  the  uncertainty  in  pH  effects  on 
absorption  are  small  and  dominated  by  the  uncertainty  due  to  ice  effects.  For  frequencies 
around  1kHz,  the  effects  could  possibly  be  significant  in  ice-free  environments  where  the 
variability  in  transmission  loss  due  to  the  uncertainty  in  SSP  is  very  small.  Study  of  the 
effects  of  pH  in  realistic  ice-free  open-ocean  scenarios  could  be  beneficial  if  the 
uncertainty  in  transmission  loss  due  to  other  factors  is  sufficiently  low. 

Since  the  pH  effects  have  been  shown  to  be  small,  and  large-scale  pH  changes  are 
projected  to  take  place  over  decades,  it  would  be  more  feasible,  but  nearly  as  beneficial, 
to  incorporate  data  from  only  periodic  pH  profiling  into  transmission  loss  models. 
Because  of  the  slow  process  of  air-ocean  gas  exchange  and  the  low  vertical  diffusivity, 
infrequent  monitoring  is  possible,  but  due  to  the  large  differences  in  summer  and  winter 
surface  boundary  conditions,  determination  of  the  structure  of  the  pH  profile  on  a 
seasonal  basis  is  recommended. 

An  unresolved  issue  is  the  pH  scale  in  use  by  different  groups.  Since  there  is  not  a 
standard  pH  scale  between  acousticians  and  chemical  oceanographers,  there  can  be  an 
offset  between  the  pH  in  the  absorption  equations  and  the  pH  measured  by  chemists,  even 
though  this  offset  would  not  affect  the  structure  of  the  pH  profile. 
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Because  of  the  issue  with  the  temporal  drift  of  the  pH  sensor,  the  manufacturer  is 
assessing  the  sensor  in  order  to  detennine  possible  corrections  to  the  measured  pH.  Once 
this  is  completed,  the  model  scenarios  could  be  repeated  and  reevaluated. 

2.  Ice  Scatter-Absorption 

The  magnitude  and  variability  of  transmission  loss  due  to  ice  scatter-absorption 
were  the  most  significant  factors  affecting  the  cumulative  TL  for  the  rays.  This  data  set 
was  from  1988,  and  the  associated  draft  values  are  larger  than  would  be  representative  for 
the  Arctic  Ocean  today.  As  sea  ice  thins  and  the  extent  declines,  it  will  be  important  to 
update  acoustic  propagation  models  and  TL  models  in  order  to  predict  accurately  how 
sound  interacts  with  the  surface/ice  boundary.  Along  with  this  updating  of  models,  it  is 
important  to  verify  them  with  empirical  data.  Thinner  ice  forms  leads  and  ridges  more 
easily,  but  the  distribution  of  keel  draft  and  spatial  frequency  can  be  expected  to  change 
from  the  thicker  ice  cover  seen  before  the  last  10  years.  Continued  efforts  to  monitor  the 
extent,  thickness  and  ice-keel  draft  of  the  Arctic  sea  ice,  and  to  adapt  predictive  models  to 
the  changing  environment  are  recommended. 

3.  Investigate  Submarine  Detection  Ranges 

To  investigate  the  effects  of  the  variability  in  T,  S,  pH,  and  ice  roughness  on 
submarine  detection  ranges,  the  model  could  be  set  up  for  given  source  noise  levels, 
source  and  receiver  depths,  environmental  noise,  and  receiver  characteristics.  Instead  of 
calculating  how  much  TL  accumulates  over  a  given  range,  and  the  factors  that  influence 
the  different  contributors,  the  model  could  determine  the  range  limit  at  which  a  given  ray 
would  no  longer  be  detectable.  To  increase  the  accuracy  and  applicability  of  this 
modeling  effort,  a  model  should  be  implemented  that  allows  an  accurate  detennination  of 
total  TL  instead  of  an  estimation  along  individual  ray  paths.  Then,  for  a  given  received 
sound  pressure  level,  and  a  known  source  level,  an  estimate  of  range  to  the  contact  could 
be  improved  by  including  the  pH  structure. 
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4. 


Model  Validation 


At  APLIS,  immediately  following  the  radial  survey,  a  group  from  the  Woods 
Hole  Oceanographic  Institute  conducted  acoustic  transmission  loss  measurements,  with  a 
source  located  at  APLIS,  and  receivers  suspended  in  the  water  column  at  each  of  the 
radial  locations  (every  9.4km).  Further  study  and  collaboration  would  be  effective  to 
compare  the  actual  TL  measurements  with  the  predicted  levels  from  the  model. 
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